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ABSTRACT 

We present an improved estimate of the occurrence rate of small planets orbiting small stars by 
searching the full four-year Kepler data set for transiting planets using our own planet detection 
pipeline and conducting transit injection and recovery simulations to empirically measure the search 
completeness of our pipeline. We identified 156 planet candidates, including one object that was 
not previously identified as a AepZer Object of Interest. We inspected all publicly available follow-up 
images, observing notes, and centroid analyses, and corrected for the likelihood of false positives. We 
evaluated the sensitivity of our detection pipeline on a star-by-star basis by injecting 2000 transit sig¬ 
nals into the light curve of each target star. For periods shorter than 50 days, we find 0.56('(g Q5 Earth- 
size planets (1 — 1.5 R©) and 0.46lQ;g5 super-Earths (1.5 — 2R0) per M dwarf. In total, we estimate 
a cumulative planet occurrence rate of 2.5 ± 0.2 planets per M dwarf with radii 1 — 4R0 and periods 
shorter than 200 days. Within a conservatively defined habitable zone based on the moist greenhouse 
inner limit and maximum greenhouse outer limit, we estimate an occurrence rate of O.lOfto'oT Earth- 
size planets and 0.12 ^q|!) 5 super-Earths per M dwarf habitable zone. Adopting the broader insolation 
boundaries of the recent Venus and early Mars limits yields a higher estimate of 0.24('(Q j!)g Earth-size 
planets and 0.21('(Q jg super-Earths per M dwarf habitable zone. This suggests that the nearest po¬ 
tentially habitable non-transiting and transiting Earth-size planets are 2.6 ± 0.4 pc and lO.Olij^ g pc 
away, respectively. If we include super-Earths, these distances diminish to 2.1 ± 0.2 pc and PC- 

Subject headings: catalogs - methods: data analysis - planetary systems - stars: low-mass - surveys 
- techniques: photometric 


1. INTRODUCTION 

In this paper, we focus on the population of planets 
orbiting small stars. Such planets are important for 
constraining the galactic census of exoplanets because 
the majority of stars in the galaxy ar e low-mass stars 
(| Henry et al.ll2006t [Winters et al.l 120151) . In addition to 
understanding the overall occurrence rate of planets or¬ 
biting low-mass stars, we would like to know how planet 
occurrence depends on factors such as planet radius, or¬ 
bital period, stellar insolation, and host star properties. 
Furthermore, small stars afford the best near-future op¬ 
portunities for detailed charact erization studies of small 
plane ts and their atmospheres (jCharbonneau fc Demin^ 
I2007D . In order to prepare for these observations, we 
would like to know the likely distance to the closest such 
targets. 

Computing the occurrence rate of small planets around 
small stars is complicated by the fact that the parameters 
of low-mass stars are more difficult to measure than the 
parameters of Sun-like stars. The main emphasis of the 
Kepler mission was the detection of planets around Sun¬ 
like stars, so the assumptions made in the construction 
of the Kepler Input Catalog (KIC) were t ailored to be 
appro priate for Sun-like stars. Accordingly. iBrown et al.l 
(1201 ID cautioned against relying on KIC classifications 
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for stars cooler than 3750K. 

Several previous studies have attempted to improve 
the KIC parameters fo r the coolest Aep/er target stars . 
In our previous paper ([Dressing Charbonneaull20I3D . 
we used the photometry provided in the KIC to 
reclassify all stars cooler than 4000K using models 
(IDotter et al.l 120081: iFeiden et al.l I2011D and a,ssump- 
tions more appropriate for low-mass stars. iGaidosI 
(I2OI3D conducted a similar analysis for the population 
of planet candidate host stars. Other authors further 
constrained the properties of particular low-mass stars 
and associated planet candidates by acquiring follow¬ 
up spec troscopic and high r esolution imaging obse r¬ 
vations (iJohnso n et al.l l2 0I2t iMuirhead et al.l ^I2allbl: 
Ballyd et al.l ~13l : iMuirhead et al.l 120131 : iSwift et all 

2013D . Recognizing the importance of characterizing the 
full target sample as well as the planet host stars in or¬ 
der to constrain the planet occurrence rate, IMann et al.l 
(I2012D acquired spectra for a subset of non-planet host 
stars. They found that the majority of bright {Kp < 14) 
Kepler target stars are giant stars (several hundred of 
which were classified as dwarfs in the KIC) but that 93% 
of fainter stars are correctly classified as dwarfs. 

In this paper, we combine the current best estimates 
of the properties of small Kepler target stars in or¬ 
der to estimate the frequency of small planets around 
small stars. Our analysis was preceded by several 
studies of the planet occurrence rate based on Ke¬ 
pler data and we adopt some techniques from the ear- 
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lier studies. In particular, we draw ut 

)on the frame- 

work established bv [Howard et al. 

([2012 

). [Fressin et al. 

(201,31). [Dressing & CharbonneaiJ ( 

2 Qli, 

[Petigura et al. 

(2013bD. and [Petigura et al.[ ( 

201,3a[). 


Working wi t h the first three quarters of Kepler data, 
[Howard et alJ l)2012[ i estimated the frequency of plan¬ 
ets around main-sequence GK stars. They found that 
the occurrence rate of planets increased sharply with de¬ 
creasing planet size and moderately with increasing or¬ 
bital period. They also found evidence for a cutoff pe¬ 
riod below which the planet occurrence rate falls off more 
quickly with decreasing period. The position of the cut¬ 
off period appeared to move outward from near 2 days 
for larger planets {Rp > SR©) to roughly 7 days for 
2 — 4 Rm planet s. 

lYoudinl (l 2011 fl used the search completeness estimates 
from [Howard et ahl (j2012fl to model the planet occur¬ 
rence rate around Sun-like stars by a joint powerlaw in 
orbital period and planet radius. He found an occurrence 
rate of 0.19 planets per star with periods shorter than 
50 days and radii larger than 2 R 0 . He also extrapolated 
outward to predict an occurrence rate of roughly three 
Earth-like planets per star with periods shorter than a 
year. 

Due to the low false positiv e rate expected for the Ke ¬ 
ple r planet candidate sa mple ([Morton fc Johnsonl[ 201 lT) . 
the [Howard et all l)2012f) and [Youdin[ ([201 ll) analyses as¬ 
sumed that all of the candidates were bona fide transiting 
planets. They further assumed that Kepler would have 
been able to detect all transiting planets with cumulative 
SNR above a set threshold of lOcr. The first assumption 
biased their occurrence rate estimates toward higher val¬ 
ues while the second assumption would have resulted in 
an underestimate if the actual search completeness were 
lo wer. _ 

[Fressin et al.[ (j2013[ l conducted a follow-up study of the 
Kepler planet occurrence rate incorporating both con¬ 
tamination from false positives and a more sophisticated 
model of pipeline sensitivity. In particular, they used a 
hierarchical approach in which they first estimated the 
population of Jupiter-size planet candidates that might 
be astrophysical false positives. They then iteratively 
determined the occurrence rate of small planets by mod¬ 
eling the fraction of larger planet candidates that might 
be masqueradi ng as smaller planet candidates in diluted 
transit events. [Fressin et al.[ (j2013[ l found a global false 
positive rate of 9.4 ± 0.9% and noted that considering 
false positives is particular important when calculating 
the occurrence rate of giant planets (6 — 22 Rq, FP rate 
= 17.7%±2.9) and Earth-size planets (0.8—1.25R 0 , FP 
r ate = 12.3% ± 3.0). 

[Fressin et al.[ ()201.‘1[ 1 also used their hierarchical model 
to estimate the completeness threshold of the Ke¬ 
pler pipeline. They found that a linear ramp model in 
which 0% of signals with SNR < 6 and 100% of signals 
with SNR > 16 were detected provided a better fit to 
the observed planet candidate population than an abrupt 
step function. Accounting for both a non-z ero false pos¬ 
itive r ate and a ramp sensitivity model, [Fressin et al.l 
(1201,81) estimated that 14.9 ± 2.4% of FGK stars host 
an Earth-size planet (0.8 — 1.25 R 0 ) with a period be- 
tween 0.8 and 50 da ys. In a similar study of Kepler data, 
[Dong fc Zhu[ ()2013[) found that roughly 20% of main se¬ 
quence stars with 5000K < Tes < 6500K host 1 — 2 R 0 


planets in period s less than 50 days. 

More recently, [Petigura et al.[ (|2013bD developed their 
own planet search pipeline in order to search for ad¬ 
ditional planet candidates around Kepler stars. Their 
TERRA pipeline uses a custom light curve detrend- 
ing algorithm based on principal component analysis 
([Petigura fc Marc^ [201211 . After search ing for planets 
around 42,557 relatively quiet GK stars, [Petigura et al.l 
([2013all found that 7.7 ± 1.3% of GK stars host small 
planets (1 — 2 R 0 ) in periods between 25 and 50 days. 
They also extrapolated to predict that 22% ± 8 % of GK 
stars host 1 — 2 R 0 planets receiving between 1 /4 and 
4 times the insolation received by the Earth. Their cal¬ 
culation incorporated a 10 % correction for false positives. 
As a be nefit of writing their own pipeline. [Petigura et al.l 
(j2013all were able to explicitly measure the completeness 
of their planet sample by injecting and attempting to 
recover transiting planets^ _ 

In a follow-up study, [Foreman-Mackev et al.l (|2014[) 
used the re ported search completa ress and planet candi¬ 
dates from fPetigura et ^ (|2013aD to rederive the planet 


occurrence rate using a hierarc hical Bayesian model. The 
ForeniamMackCTet_aH (|2014[ ) analysis differed from the 


Petigura et al.l (2013aD analys is in two key aspects: (1) 


Foreman-Mackev et al.[ (|2014f) considered measurement 


errors in the stellar and transit parameters and ( 2 ) they 
did not assume that the planet occurrence rate was flat 
in log period, instead using a flexible Gaussian process 
to model the occurrence rate assuming a smooth func - 
tional form. As a result, [Foreman-Mackev et al.[ (j20l3) 
found an occurrence rate of potentially h abitable Earth- 
size pla nets three times lower than the [Petigura et al.l 
([2013al) estimat e. _ 

[Silburt et al.l (|201,5[ ) considered a sample of 76,711 Ke¬ 
pler target stars with radii of 0.8 — 1.2 R 0 and es¬ 
timated the search completeness using the reported 
Combined Differential P hotometric Precision (CDPP, 
[Christiansen et ahl [2012f) . They employed an iterative 
simulation to investigate the dependence of the planet 
occurrence rate on planet radius without subdividing 
the data into bins and accounted for errors in the 
planet radii. For p lanets with periods of 20-200 days, 
[Silburt et al.[ (|2015[) reported that the occurrence rate 
is higher for planets with radii of 2 — 2.8 R 0 than for 
smaller or larger planets. In total, they estimated that 
a typical Sun-like star hosts 0.46 ± 0.03 planets with pe¬ 
riods of 20 -200 days and radii o f 1 — 4Rm. In agree¬ 
ment with [Petigura et al.[ (j2013aD . [Silburt et al.[ (j201,5f ) 
noted that the planet occurrence rate is flat in log pe¬ 
riod. Within a broad hab i table zone extending from 
0.99-1.7 AU, [Silburt et al.1 (|2015D estimated an occur- 
rence rate of 0.064to,oii small (1 — 2 R 0 ) planets per 
Sun-like star. 

F ocusing specifically on Kepler ’s smallest target stars, 
we (jDressing fc Charbonneadl201,3[ ) estimated an occur¬ 
rence rate of O.9Oto03 planets per star for 0.5 — 4 R 0 
planets with periods shorter than 50 days. We based 
our previous analysis on Q1-Q6 Kepler planet candidate 
list and assumed that Kepler detected all planets with 
cumulative SNR > 7.1g. Using cons ervative habitable 
zone limits from [Kasting et HI (jl99.3f ). we estimated an 
occurrence rate of O.lbj’jg jg potentially habitable Earth- 
size (0.5 — 1.4 R 0 ) planets per small star. [Kopparar^ 
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(|2013ll then revised this estimate to 0.481 q' 24 planets per 
star using the broader updated hab itable zone bound¬ 
aries from iKopuarapu et alJ (l2013bf) . His result agreed 
well with the estimate of 0.46to'i5 potentially habitable 
0.8 — 2Rm planets p er star from IGaidosI (|2013fl . Un¬ 
like [KopE^ra^ (j2013[ l. [Gaido^ (j2013l l adopted habitable 
zone bou ndaries correspond ing to the 50% cloud cover 
c ase fromiSelsis et a.1.1 (120071) . 

iMorton fc SwiftI (1201411 adopted a slightly different 
technique to estimate the frequency of small planets 
around small stars. They assumed that the planet 
radius distribution is independent of orbital period 
and modeled each planet using a weighted kernel den¬ 
sity estimator when computing the occurrence rate. 
They found that the oc c urrenc e rate estimat e s from 
iDressing fc Charbonneaul (|2013f) . iKopoara^ (1201311 . 
and IGaidosI (|2013l l for planets smaller than 1.4 R® should 
be increased by an additio nal incompletenes s facto r of 1.6 
if the assumption made bv IMorton fc Swif3 (|2014ll about 
the period-independence of the planet radius distribution 
is correct. _ 

Like lSilburt et al.l (|2015ll . iGaidos et ahl (1201411 used an 
iterative simulation to estimate the planet occurrence 
rate, but they elected to focus on stars cooler than 
4200K. For or bital periods of 1 — 180 days and radii 
of 0.5 — 6Re, IGaidos et al.l (|2014ll calculated a cumu- 
lative occurrence ra te of 2.01 ±0.36 planets per M dwarf. 
IGaidos et al.l (|2014ll also remarked that the planet oc¬ 
currence rate is highest for planets with radii of approx¬ 
imately 1 Re and lower for larger and smaller planets. 

The frequency of potentially habitable planets around 
small stars has also been estimated from radial veloc¬ 
ity surveys. Based on six years of observa tions with the 
HARPS spectrograph. iBonfils et all ()2013ll estimated an 
occurrence rate of O.dl^Q ^g potentially habitable plan¬ 
ets per M dwarf. Their definition of “potentially hab¬ 
itable” encompassed planets with 1 < Msini < 10 
within the “early Mars” and “recent Ven us” bound¬ 
aries of the habitable zone p resented in iSel sis et all 
(|2007ll . A subsequent study bv iRobertson et al.l ( 2014f) 
reveal ed that GJ 581d. o ne of the two planets upon 
which IBonfils et al.l (j2013ll based their occurrence rate 
estimate, is like l y a m anifestation of stellar activity. 
iRobertson et al.1 (j2014H reported a revised occurrence 
rate of 0.33 potentially habitable planets per M dwarf. 
The updated RV-based estimate is more similar to the 
estimates based on Kepler data, but accurately deter¬ 
mining an o ccurrence rate with only a single planet 
(GJ ( 167Gc. lA nglada- Escude et al.l 120121 : IBonfils et al.l 
l2013t iDelfosse et al.11201311 is challenging. Additionally, 
direct comparison of planet occurrence estimates from 
RV and transit surveys is complicated by the need to 
employ a compositional model to translate planet masses 
into radii. 

In this paper, we implement the following improve¬ 
ments to refine our 2013 estimate of the frequency of 
small planets around small stars: 


We use the full Q0-Q17 Kepler data set. 


• We utilize archival spectroscopic and photometric 
observations to refine the stellar sample. 


• We inspect follow-up observations of planet host 
stars to properly account for transit depth dilution 
due to light from nearby stars. 

• We apply a correction for false positives in the 
planet candidate sample. 

• We incorporate a more sophisticated treatment of 
the habitable zone. 

In Section we describe the selection of our stellar sam¬ 
ple, which includes some stars whose parameters have 
been characterized spectroscopically via follow-up obser¬ 
vations. We explain our planet detection pipeline in Sec¬ 
tion [3] and our procedure for vetting candidates in Sec- 
tiondl We present light curve fits for the accepted planet 
candidates in Section [51 In Section [51 we assess the com¬ 
pleteness of our pipeline. We then estimate and discuss 
the planet occurrence rate in Section [7] before concluding 
in Section [51 


2. STELLAR SAMPLE SELECTION 


We selected our stellar sample by first downloading a 
table of all 4915 stars with TeS < 4000 and logg > 3 
from the Ql-16 Kepler SteWai Cat alog on the NASA 
Exoplanet ArchivJ^ (lAkeson et al.11201311 . This cata¬ 
log is described in iHuber et al.l (|2014D and combines 
the best estimates available for each star from a va¬ 
riety of photometric, spectroscopic, and asteroseismic 
analyses. The properties for the stars in the down¬ 
loaded sa mple were primar il y determined from pho¬ 
tome t ry (iBrown et al.l I20ITI: IDressing fc Charbonneaul 
I 2 OI 3 I : iGaidod l2013t iHuber et al.l l2014ll ! but 2% of 


the sample had s pectroscopically-derived parameters 
(iMann et al.l 120121: Mui rhead et all l2012al: iMann et al.1 
I^OldbUMartfn et al. I2013h . For the majority of the stars 


in the sample (79%), t he stellar parameters were draw n 
from our 2013 analysis ([Dressing fc Charbonneaiill2013[l . 

Some of the stars in the downloaded subset had light 
curves indicative of binary stars, variable stars, or en¬ 
hanced spot activity. In addition, some of the stars were 
observed only for a small number of days. In order to 
accurately estimate the planet occurrence rate for small 
stars, we wanted to select the subset of stars with the 
highest search completeness. We therefore performed the 
following series of cuts on the sample. 

First, we counted the number of timestamps for which 
each star had “good” data (i.e., not flagged). We re¬ 
jected all 2101 stars with fewer than 48940 unflagged 
long cadence data points. Since Kepler obtained long 
cadence data using 29.4 minute integration times, this 
cut requires 1000 days of data. One of the main goals 
of this paper is to measure the occurrence rate of poten¬ 
tially habitable planets and we wanted to ensure that Ke¬ 
pler would have been able to observe multiple transits of 
planets within the habitable zones (HZ; see Section [77311 
of the stars in our final sample. For refe rence, the median 
orbital period at the outer edge of the iKopparapu et al.1 
(|2013bll HZ for the stars in our final sample is 131 days 
and the longest period at the outer HZ is 207 days. 

Seco nd, we removed 63 stars that iMcQuillan et al.1 
(|2013[) categorized as likely giants based on their stochas¬ 
tic photometric variability. The affected stars have red 


• We explicitly measure the pipeline completeness. 


® http://exoplanetarchive.ipac.caltech.edu 
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Fig. 1.— Radii and stellar effective temperatures of the stars in 
our final selected subsample (blue) compared to the full sample 
initially downloaded from the NASA Exoplanet Archive (black). 
The red lin e is an empirical rel ation betwee n eff ective temperature 
and radius IIMann et al.ll2013al . see Section l7.4t . 

colors (median J — H = 0.83) consistent with their re¬ 
vised classification as giants. Although two of the stars 
had re vised classifications from iDressing &: Charbonneaul 
(|2013f) . the remaining 61 had parameters from the Kepler 
Input Catalog. 

For reference, we checked whether any of our target 
stars were known eclipsing binaries by consulting the Ke¬ 
pler Eclipsing Binar y Catalog. We exam ined both the 
published Version 2 dSlawson et ahll^OllD and the online 
beta version of the Third Revision^ Eleven of the target 
stars were listed in both versions of the catalog, six stars 
were listed in Version 2 only and three stars were listed 
in Version 3 only. 

Six of the twenty targets with matches in the Eclips¬ 
ing Binary Catalogs were listed in the NASA Exo¬ 
planet Archive as false positive systems (KID 5820218 
= KOI 1048, KID 6620003 = KOI 1225, KID 8823426 
= KOI 1259, KID 9761199 = KOI 1459, KID 9772531 
= KOI 950, and KID 10002261 = KOI 959). Two were 
listed as planet candidate host stars (KID 5384713 = 
KOI 3444 and KID 11853130 = KOI 326 3). Confirmed 
giant planet KOI 254.01 (KID 5794240, l.lohnson et al.l 
I20I2D was also included as an EB match because of the 
very large transit depth. As evidenced by the presence 
of KOI 254.01 in the EB catalogs, the catalogs contain 
both actual EBs and likely planets. Accordingly, we did 
not remove the twenty targets with matches in the EB 
catalogs from our target sample. 

We then detrende d al l of the light curves as described 
in detail in Section [ST] using smoothing lengths of 500, 
1000, and 2000 minutes. We constructed a histogram 
of the flux distributions for each of the detrended light 
curves and measured the of a fit to a Gaussian flux 
distribution. The flux distribution of a well-behaved sin¬ 
gle star should be Gaussian after detrending, but the 
flux distribution of an eclipsing binary can appear bi- 
modal. We therefore flagged for visual inspection all 
353 stars for which reduced > 3 for any of the de- 

keplerebs.villanova.edu 


trended light curves. We also measured the standard 
deviations ctsooi ctioooj o '2000 of the three detrended light 
curves for each star and took the ratios of the standard 
deviations of light curves detrended using different me¬ 
dian filters. We flagged 42 stars for which any of the 
ratios crsQo/criooo, 0 - 500 /cr 2 ooo, or crioQo/o- 2 ooo were below 
0.8. This cut was designed to pick out light curves for 
which the detrending algorithm failed to remove longer 
timescale variability. 

Finally, we visually inspected th e detre n ded and PDC- 
MAP photometry ([Smith et al.l I2012t IStumne et al.l 
I2012D for the 395 flagged stars with > 1000 days of data 
and large or standard deviation ratios. We rejected 

207 stars with highly variable detrended light curves or 
strong indications of classification as an eclipsing binary. 
We also rejected KID 12207013 because of the unusual 
light curve morphology displayed in certain quarters. We 
accepted the remaining 187 flagged stars, increasing our 
final selected sample to 2543 stars. We note that the 

208 stars that were rejected during the visual inspection 
stage are unlikely to harbor detectable planet candidates 
exactly because their light curves are highly variable. 
Similarly, our injection tests would likely recover only 
a small fraction of any planets injected into their light 
curves. Accordingly, the exclusion of these stars from 
our stellar sample has negligible effect on our estimated 
rates of planet occurrence. 

Our final selected sample of 2543 stars is compared to 
the initial downloaded sample in Figure |TJ The temper¬ 
ature range of the sample extends from 2661K to 3999K, 
with a median stellar effective temperature of 3746K. 
The median stellar radius is 0.47 R 0 and the stars have 
radii spanning from 0.10 R© to 0.64 R©. The metallicity 
range is [Fe/H] = -2.5 to [Fe/H] = 0.56, with a slightly 
sub-solar median metallicity of [Fe/H] = —0.1. How¬ 
ever, most o f the metallicity estimates were d erived from 
photometry ([Dressing fc Charbonneaull2013l) and are not 
well-constrained. The brightest star in the sample has a 
Aep(er magnitude Kp = 10.07, but the median bright¬ 
ness is Kp = 15.5. The faintest star has Kp = 16.3. The 
sample contains 100 known planet (candidate) host stars 
with 83 planet candidates and 80 confirmed planets. 

3. PLANET DETECTION PIPELINE 

The first step in our planet detection pipeline was to 
clean the light curves to prepare them for the transit 
search. Next, we searched each light curve sequentially 
for planets, allowing the code to detect multiple plan¬ 
ets per star when warranted by the data. We then 
accepted or rejected putative detections using the vet¬ 
ting procedure described in Section Uj and checked that 
the automatically accepted transits were not ephemeris 
matches with other KOIs. Finally, we visually inspected 
all surviving candidates and reviewed all available follow¬ 
up analyses produced by the Kepler team and the com¬ 
munity. We discuss each step of the planet detection 
pipeline in more detail in the following sections. 

3.1. Preparing the light curves 

We obtained all available long cadence data for each 
target via anonymous ftp from the MAST0 We then 

® http://archive.stsci.edu/kepler/publiclightcurves.html 
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Fig. 2.— Illustration of the detrending process using a section 
of the light curve of KID 5531953 (KOI 1681). Top: PDC-MAP 
flux versus time. Middle: Normalized flux versus time. Bottom: 
Flux detrended using a 2000 minute filter versus time. 

excluded all data points flagged as low quality and de¬ 
trended each quarter of data independently. We pro¬ 
duced three detrended versions of each light curve using 
a running sigma-clipped mean filter with widths of 500, 
1000, or 2000 minutes. When calculating the mean, we 
excluded all points more than Scr away from the median 
value of the light curve within the filtered region. We 
then divided the flux data by the smoothed light curve 
to obtain a detrended, normalized light curve for that 
quarter. Figure [5] provides an illustration of our light 
curve detrending process. 

Next, we searched for data gaps and anomalies within 
the detrended light curves. We defined a data gap as 
> 0.75 days of missing photometry. Data gaps are fre¬ 
quently accompanied by sharp increases or decreases in 
flux that can confound searches for planets. We removed 
these events by excising all data points within 1 day of 
the start of a data gap or less than 3 days after the end 
of a data gap. 

3.2. Searching for Transiting Planets 

We designed our planet search algorithm to take ad¬ 
vantage of the known stellar properties for each target 
star. First, we predicted the expected transit duration 
as a function of orbital period base d on the mass and 
radius of the target star (jWinnll^OKIf) . We initially com¬ 
puted the transit duration for a central transit of a planet 
in a circular orbit, but we reduced the minimum dura¬ 
tion considered by a factor of 15 to account for graz¬ 
ing transits and eccentric orbits. We then constructed a 
box-fitting least squares (BLS) periodogram for each of 
twelve logarithmically spaced intervals between 0.5 and 
200 days. Our main planet search program was writ¬ 
ten in IDL, but the BLS portion was implemented using 
the Fortran package and IDL wrapper provided by Scott 
Fleming 0 

We used different boundaries for the transit “duty cy¬ 
cle” (the ratio of the transit duration to the orbital pe¬ 
riod) for each period range based on our predicted tran- 

® http: //www .personal .psu. edu/users/s/w/swf 13/SGE/clio .html 


sit durations. For each search, we used a light curve de¬ 
trended with a smoothing filter of 500, 1000, or 2000 min¬ 
utes. The choice of light curve was set by the expected 
transit duration. For predicted transit durations shorter 
than 200 minutes, we selected the shortest smoothing fil¬ 
ter such that the expected transit duration was less than 
one tenth of the smoothing window. For transit dura¬ 
tions longer than 200 minutes we used the 2000 minute 
filter. 

We then determined the signal detection efficiency 
(SDE) for each possible sign al in the cornposite peri¬ 
odogram using Equation 6 in iKovacs et al.l (|2002[ ). We 
checked whether any peaks had SDE > 6 and stopped 
searching if no peaks were above the threshold. If peaks 
were detected, we ranked the peaks in order of decreasing 
SDE. 

Starting with the most significant peak, we re-ran the 
BLS algorithm considering only periods close to the pe¬ 
riod of the identified peak. In these high-resolution runs, 
we considered transit durations (rdur) as short as l/30th 
of the expected transit duration of a planet in a circu¬ 
lar orbit to increase the chance that our code would be 
able to recover planets in grazing or eccentric orbits. We 
used these higher resolution BLS periodograms to deter¬ 
mine the epochs of the putative transit events. We then 
ran a Mon te Carlo analysis to fi nd the preliminary tran¬ 
sit model (|Mandel &: Agoll[2002li that best described the 
candidate event. 

In our Monte Carlo analysis, we allowed the transit 
center to shift by 2 hours (up to a maximum of I/200th 
of the orbital period for short-period events). For each 
choice of transit center we generated a new version of the 
detrended light curve by dividing the raw PDC-MAP 
Kepler photometry by a straight line fit to the pho¬ 
tometry immediately preceding and following each pu¬ 
tative transit. Specifically, we considered data points 
more than one and less than 3.5 expected full transit 
durations away from the putative transit center. We as¬ 
sumed circular orbits and estimated quadratic limb dark¬ 
ening parameters from the Tgg and log g of the target star 
by interpolating b e tween the coefficients determined by 
iClaret fc BloemenI (|20lH ). We considered a/i?* between 
50% and 200% of the expected value for the trigger or¬ 
bital period, impact parameters between 0 and I, and 
Rp/Rt: as large as the square root of the depth of the 
trigger event. 

In some cases, the highest peak in the periodogram 
was actually at a harmonic of the true planet period, so 
we repeated the Monte Carlo transit fitting analysis at 
i and n times the trigger period for n = 2 — 7. We then 
selected the period for which the compared to a 

straight line fit was maximized. We rejected all putative 
transit events for which none of the models were preferred 
at 5cr and recorded the parameters of the best-fit model 
in other cases. We further refined the planet parameters 
during the vetting stage as described in Section [H 

We then repeated the process described in the previous 
three paragraphs to fit the next highest peak in the peri¬ 
odogram. If the preceding transit fit had been accepted, 
then we excised the data near transit prior to fitting the 
next peak. When all peaks above the threshold level 
were exhausted, we generated a new periodogram using 
only the out-of-transit data and reran the peak identi- 
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fication and transit model fitting process with the new 
periodogram. The code automatically stopped searching 
for planets when no peaks with SDE > 6 were found, 
when none of the transit models for the identified peaks 
were accepted, or when the code had completed three it¬ 
erations of searching for planets. (Note that multiple 
planets could be detected in a single round of searching.) 

4. VETTING 

Our transit detection pipeline identified 3111 puta¬ 
tive transit events associated with 534 stars. Some of 
those signals might have been systematics or astrophys- 
ical false positives instead of bona fide transiting planet 
candidates. Accordingly, we performed a series of cuts 
to select the events consistent with transiting planets. 
First, we visually inspected the candidate transit events 
to identify signals that were not clearly associated with 
spacecraft systematics or stellar activity. Of the 3111 
candidate events, 511 events survived initial visual in¬ 
spection. The 2600 events rejected at the visual inspec¬ 
tion stage displayed morphologies consistent with classi¬ 
fication as spacecraft systematics or sinusoidal brightness 
variations indicative of starspots rather than transiting 
planets. The 511 signals surviving visual inspection were 
associated with 246 unique host stars. 

Some of the candidate events were harmonics of signals 
detected at integer multiples of the true period. We then 
ranked the accepted signals for each host star in order of 
decreasing as calculated during the detection phase 
and iteratively fit and excised the transits of each signal 
in order to ensure that the lower Ay^ signals were not 
simply harmonics of the strongest signals. We rejected 
all signals with resulting Ay^ below 5 ct. This “sequential 
vetting” step reduced the number of candidate events to 
323 possible transits for 246 unique stars. 

Next, we conducted a second, more intensive round of 
visual examination for the remaining candidate events. 
We compared the shapes and depths of odd and 
even transits, checked for the appearance of secondary 
eclipses, considered the depth of the putative transit rel¬ 
ative to other possible features at the same orbital pe¬ 
riod, and investigated whether the putative transit events 
were dominated by a small number of deep events. Af¬ 
ter visually vetting the candidates, we checked whether 
putative signals had been previously identified as false 
positives and excluded all such events. We rejected 180 
signals (including 9 known false positives) and accepted 
143 signals associated with 97 unique stars. 

During the vetting stage, we noticed that the phase- 
folded light curve for KOI 2283.01 exhibited two transit¬ 
like events with markedly different depths when folded 
to the 17.402 day orbital period listed in the Ql-16 KOI 
catalog. We therefore rejected KOI 2283.01 as a blend 
containing an eclipsing binary. This interpretation is 
consistent with the large observed centroid source offset 
shift of 5.3 ct. 

We then performed a final search for additional planets 
in the systems in which a previously detected planet had 
survived the vetting process. We executed this search 
by phase-folding the detrended light curve on the or¬ 
bital periods of all accepted signals and removing all data 
points within plus or minus one best-fit duration of tran¬ 
sit center prior to re-running the search process detailed 
in Section 13.21 In all cases we used the light curve that 


TABLE 1 

New Candidate Accepted By Our Pipeline 


KID 

KOI 

P 

(days) 

to 

(days) 

Rp 

(R-©) 

a 

(AU) 

Ax2 

5531953 

1681 

21.913843 

17.036402 

1.03 

0.16 

41.4 


had been detrended using a 2000 minute filter. The mo¬ 
tivation for repeating this search after the first round of 
vetting was that uncertainties in the initial periods, du¬ 
rations, and transit centers of the accepted signals might 
have limited the effectiveness of the clipping performed 
in the initial search. 

Our second round search revealed 104 candidate signals 
for 33 stars. We vetted the signals using the same vetting 
pipeline as in the first round search and accepted 15 ad¬ 
ditional candidate transiting planets associated with 12 
stars. Next, we excised the transits of the signals ac¬ 
cepted in the second round and performed a third round 
of transit searches. No additional signals were accepted 
during the third round. The full sample of 156 planet 
candidates included 143 signals associated with 97 stars 
revealed in the first round of searching and 15 signals 
associated with 12 stars revealed in the second round. 

We compared the periods, P, and epochs, to, of the ac¬ 
cepted planet candidates to the catalog of known eclips¬ 
ing bmurieSj_geriodic_raTiable stars, and KOIs compiled 
by iCoiighlin et al.l (j2014[l . We excluded the host star 
from the match process in order to avoid matching a sig¬ 
nal to itself. We did not find any corresponding signals 
within our specified match tolerances of 

l^inatch - < min (2 hr, 0.001 x P) 

and (1) 

|to, match -to\< min (4 hr, 0.001 x P) 

4.1. New Planet Candidate 

The majority of the accepted planet candidates corre¬ 
sponded to signals previously identified as KOIs. We 
found that 155 putative planet candidates had peri¬ 
ods and epochs matching t hose of known planet cand i¬ 
dates or confirmed planets (iBonicki et al.l[2?)Tol l2011al lbl: 
iBatalha et al.ll201^ iBurke et al.ll2014[ l. We accepted one 
new signal in a system with previously known KOIs. 

The KOI 1681 (KID 5531953) system contains three 
known planet candidates with periods of 6.51, 1.99, and 
3.53 days. Our pipeline detected a planet candidate in 
the system with a radius of IR 0 and a period of 21.9 
days, roughly 11 times the 1.99 day orbital period of 
KOI 1681.02. We describe additional planet properties 
in Table [TJ As shown in Figure [31 the transit signal is 
still visible in the light curve after the transits of the 
other three planets have been removed, suggesting that 
this is a new planet candidate rather than an alias of 
KOI 1681.02. The signal was later identified in the Ql- 
17 DR 24 Kepler pipeline rurQ and is now listed as planet 
candidate 1681.04. 

4.2. Accounting for Transit Depth Dilution 

^ http://tinyurl.com/kepler-dr24 


















The Occurrence Rate of Potentially Habitable Planets Orbiting M dwarfs 


7 


15 L Accepted Signals: 4 
Known KOIs: 3 

10 - 

5 = ^ 81.01 
0 


P (Days) 


1.0000 

0.9998 

0.9996 

0.9994 




3 4 

Folded Time 










1,0004 
1.0002 
1,0000 
0,9998 
0,9996 


: Even 



P:3.53 j 

0.05 0,10 0.15 0.20 0.25 

Folded Time 



Fig. 3. — Transit signals detected in the KID 5531953 system. Top left: Period and epochs of all four signals identified by our pipeline 
(circles) and the known planet candidates in the system (marked by text). Top right: Transit times for each of the four planet candidates. 
Second row from top: Detrended fiux versus time folded to the 6.9 day period of KOI 1681.01. The left panel displays the full binned 
phase-folded light curve and the right panel shows a zoomed-in view near transit center. The orange line marks the transit center. The 
light gray and dark gray lines show the binned phase-folded light curve for only the odd and even transits, respectively. We excised data 
points between the vertical gray lines before folding the data to the period of the next planet. Third row from top: Same as second row 
but for the 1.99 day period of KOI 1681.02. The transits of KOI 1681.01 are not included. Fourth row from top: Same as previous row 
but for the 3.53 day period of KOI 1681.03. The transits of KOI 1681.01 and KOI 1681.02 are not included. Bottom row: Same as second 
row but for the 21.9 day period of the new signal detected by our pipeline. The transits of KOI 1681.01, 1681.02, and 1681.03 are not 
included. 


For the 155 known KOIs in our planet candidate sam- 
pie, we inspected the DV reports prepared by the Ke¬ 
pler team and all publicly available follow-up data to 
check for signs that the KOIs were false positives. As 
described below, we learned that several of the planet 
host stars in our sample have stellar companions at sep¬ 
arations within 1", Accordingly, the measured transit 
depths for those planet candidates would have been di¬ 
luted by the additional light in the aperture. 

Two of those syste ms (KOI 1422 and KOI 2626) were 
well characterized bviCartier et al,l (|2014f) . In their anal¬ 
ysis, iCartier et all (j2014H determined stellar parameters 
for the double star system KOI 1422 and the triple star 
system KOI 2626 using HST WFC3/UVIS photometry. 
They were unable to constrain which of the stars hosted 
the associated planet candidates, but they were able to 
provide revised estimates for the radii and orbital param¬ 
eters of the associated planet candidates for each choice 
of host star. In our analysis, we therefore chose to repre¬ 
sent the planet candidates in these systems using “frac¬ 
tional planets” orbiting each of the possible host stars 
rather than assuming that the planet candidates orbit 
the system primaries or excluding them from the analy¬ 


sis. 

For the remaining systems with close stellar compan¬ 
ions we did not have sufficient information about the 
companion star to model planet candidates orbiting each 
star in the system. Instead, we corrected for the transit 
depth dilution by multiplying the estimated planet ra¬ 
dius by the correction factor c = -1/2.512“^^ -|- 1, where 
AK is the difference in Kepler magnitudes between the 
apparent magnitudes of the target star and companion 
star. In the worse case scenario of an equal-brightness 
binary, the correction increases the estimated radius of 
the planet candidate by roughly 40%. 

We applied correction factors for three systems: 
KOI 605 (41%), KOI 3010 (41%), and KOI 3284 (2%). 
In the KOI-605 system, which contains two candidates, 
D. Ciardi’s Keck/NIRC2 adaptive optics imaged re¬ 
vealed that the system consists of two stars separated 
by less than O'.'l with nearly equal brightness in the Ke¬ 
pler bandpass. 

D. Ciardi also acquired Keck/NIRC2 observationtQ of 

® https://cfop.ipac.caltech.edu/edit_obsnotes.php?id=605 

® https://cfop.ipac.caltech.edu/edit_obsnotes.php?id=3010 
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KOI 3010 showing that the system is a close binary with 
a separation of 0"3. The two stars appear to have nearly 
equal brightnesses. For KOI 3284, Keck/NIRC2 and 
Gemini/DSSI imagetO revealed a companion 3.56 mag¬ 
nitudes fainter than the target star at a separation of 
0"4. 

Our correction procedure explicitly assumed that any 
associated planet candidates orbit the target star, but 
they might actually orbit the companion star. If the 
planet candidates do indeed orbit the companion star, 
then the radii of the planet candidates will need to be 
reevaluated once t he properties of th e companion star 
are established isee lCiardi et al.ir2015l for a detailed dis¬ 
cussion). In most cases, the available photometry was 
insufficient to determine whether the nearby companion 
is physically associated with the target star or to con¬ 
strain the properties of the companion star. 

4.3. False Positive Correction 

In addition to correcting for transit depth dilution in 
systems with nearby companions, we incorporated a gen¬ 
eral false positive correction to account for the possibility 
that some of the smaller transiting planets in our sample 
might be diluted eclipses or tr ansits of larger plane ts. We 
therefore consulted Table 1 of iFressin et al.l (j2013l l to de¬ 
termine the false positive probability (FPP) for a planet 
with a given radius. We apply this radius-dependent cor¬ 
rection to the planet occurrence map derived in Section[3 

4.4. Known Planet Candidates Missed hy Our Pipeline 

Our sample of accepted candidate events included all 
but 7 of the 161 known planet candidates and confirmed 
planets meeting our sample cuts of planet radii larger 
than 0.5 Re and orbital periods shorter than 200 days. 
We list the missed candidates in Table [2] We note that 
I Swift et ahl (|2015D rejected one of the missed candidates 
(KOI 1686.01) as a possible planet because the phase 
folded light curve did not display a convincing transit 
event and that the reported disposition in the NASA 
Exoplanet Database was later changed to False Posi¬ 
tive. Three of the missed candidates are in the same 
system (KOIs 3444.01, 3444.03, and 3444.04) and none 
of the 7 missed candidates produced accepted peaks in 
the BLS periodograms. Although we were reassured that 
our pipeline recovered most of the previously detected 
planet candidates, our goal was not to reproduce the Ke¬ 
pler planet candidate list but to design a single pipeline 
that could be used to both search for planets and char¬ 
acterize pipeline completeness. Thus we do not consider 
these additional 7 KOIs in our analysis below. 

5. PLANET PROPERTIES 

Accurate planet radius estimates are a key ingredient 
in the planet occurrence calculation. We therefore re¬ 
fined the preliminary transit parameters found in Sec¬ 
tion 13.21 by conducting a Bayesian Markov Chain Monte 
Carlo (MCMC) an alysis with a Metropolis -Hastings ac¬ 
ceptance criterion (|MetroDolis et al.l 119531 1. We varied 
the orbital period P, epoch of transit center to, planet-to- 
star radius ratio semimajor axis to stellar radius 

ratio a/i?* and impact parameter b. We assumed that all 

https: //cf op. ipac. caltech.edu/edit_obsnotes .php?id=3284 


TABLE 2 

Known KOIs Missed By Our Pipeline 


KID 

KOI 

Kepler 

name 

mag 

{Kp) 

P 

(days) 

Rp 

(R©) 

Kepler 

SNR^ 

6149553 

1686.01*= 


15.89 

56.87 

1.3 

5.2== 

8890150 

2650.02 

395b 

15.99 

7.05 

1.1 

12.4 

9605552 

3102.01 


15.98 

9.32 

1.0 

6.1 

5384713 

3444.01 


13.69 

12.67 

1.0 

14.7 

5384713 

3444.03 


13.69 

2.64 

0.6 

12.2 

5384713 

3444.04 


13.69 

14.15 

0.8 

9.3 

2986833 

4875.01 


15.78 

0.91 

1.0 

10.2 


^ As reported in the Cumulative KOI table at the NASA Exo¬ 
planet Archive on 22 November 2014. 

^ As of 22 March 2015, the reported disposition of KOI 1686.01 at 
the NASA Exoplanet Archive has been changed to False Positive. 

The transit SNR for KOI 1686.01 was not reported in the Cu¬ 
mulative KOI table, the Q1-Q16 table, or the Q1-Q12 table. 
This is the value from the Ql—Q8 table. The value in the Q1-Q6 
table was 7.6. 

of the orbits were circular and fixed quadratic limb dark¬ 
ening parameters to the values predic ted from the stellar 
temp eratures and surface gravities (|Claret &: Bloemeiil 
1201111 . We conducted some of the planet fits using short 
cadence Kepler data to better constrain the shape of 
transit during ingress and egress. 

For each planet candidate, we ran N chains starting 
at initial positions set by perturbing the initial solution 
found during the detection and validation process by up 
to 5(7 in each parameter. We manually adjusted the 
step sizes for each parameter such that the acceptance 
fractions were between 10-30%. We ran each chain for 
at least 10^ steps before initiating periodic convergence 
tests by calculating the Celman-Rubin potential scale 
reduction factor R for each parameter (|Celman et al.l 
1200411 . We terminated the chains when R < 1.05 for 
all parameters and then accounted for “burn-in” by re¬ 
moving all steps taken prior to the point at which the 
likelihood first became higher than the median likelihood 
of the chain. After merging the chains, we adopted the 
median values of each parameter as the best-fit value and 
assigned errors encompassing the 68% of values nearest 
to the chosen best-fit value. We provide the best-fit pa¬ 
rameters for each detected planet candidate in Tabled 
and display them in Figured! 

Several of the planet candidates in our sample ex¬ 
hibit large transit timing variations and were poorly 
fit by the method described above. In those cases, we 
adopted the transit parameters found in previous stud¬ 
ies using combined fits to the individual transit times 
and the planet properties. For the candidates d i splay - 
ing TTVs, we incorporated fits from iRowe et al.l (|2014[ 1 
for KOIs 248.01, 248.02, 248.03, 248.04, 314.01, 314.02, 
314. 03, 448.01, 448.02, 886 .01, 886.02, and 886.03 and 
from lloannidis et al.l (j2014| l for KOIs 676.01 and 676.02. 

In addition, we adopted the light curve parame¬ 
ters from the 2 January 2015 version of the NASA 
Exoplanet Archive for 5 candidates: 961.01, 961.03, 
1681.01, 2329. 01, and 32 6 3.01. We also adopted 
the fits from IRowe et al.l (|2014[ 1 for KOIs 254.01, 
430.01, 430 . 02, 43 8.01, 775.02, 868.01, and 961.02, from 
iSwift et all (Iml) f or KO Is 1902.01 and 3444.02 and 
from iCartier et al.1 l|Mll for KOIs 1422.01, 1422.02, 
1422.03, 1422.04, 1422.05, and 2626.0 1. The KOIs with 
parameters from ICartier et al.l (I2014D were detected in 


















































The Occurrence Rate of Potentially Habitable Planets Orbiting M dwarfs 


9 




Fig. 4.— Radii of the planet candidates detected by our pipeline versus orbital period (Left) or insolation flux {Right) with Icr errors. In 
most cases, we refit the planet parameters by conducting an MCMC analysis fitting transit models to the short cadence (cr imson points) 
or lon g cadence (teal points) Kepler data.. For the remai ning 34 planet can didates, we a dopted transit par ameters from IHartier et al.l 
II2014I . purple points'). Iloannidis et al.l II2014I . brown pointsl. IRowe et al.l II2014I . black points'). ISwift et al.l II2015I . blue points), or the NASA 
Exoplanet Archive (gray points). In all cases, the errors on the planet properties incorporate uncertainties in both stellar and transit 
parameters. The green boxes indicate the boundaries within which we report planet occurrence rates and the arrows in the right panel 
mark several variations of the habitable zone as explained in the legend. The errors for KOIs 1422.01, 1422.02, 1422.03, 1422.04, 1422.05, 
and 2626.01 appear particularly large because we accounted for the possibility that any of the three stars in the KOI 2626 system or either 
of the two stars in the KOI 1422 system harbor the transiting planets (see Section 14.21 for details). KOI 254.01 has an estimated radius of 
IIR0 and therefore does not appear on these plots. 


multi-star systems in which the identities of the host stars 
are unknown. As described in Section 321 'we accounted 
for all possible system configurations by using fractional 
planets distributed around each of the possible host stars. 

As depicted in Figure IH we found that 
149 of the accepted planet candidates had re¬ 
vised radii 0.5 <i?p< 4 R 0 and orbital periods 
0.5 < P < 200 days. Of the remaining candidates, one 
had a shorter period (KOI 961.02, P = 0.45 days), 
one was too small (KOI 5692.01), and six were too 
large (KOI 254.01, 868.01, 901.01, 902.01, 1176.01, and 
3263.01). KOI 868.01 also has an orbital period longer 
than 200 days. 

6. PLANET INJECTION PIPELINE 

In order to accurately measure the planet occurrence 
rate based on the results of our planet search, we needed 
to know the completeness of our planet candidate list. 
We measured the completeness of our planet detec¬ 
tion pipeline by injecting transiting planets into the 
PDC-MAP light curves, detrending them, and running 
the detrended light curves through our detection algo¬ 
rithm. We did not introduce the signals at the pixel level 
and we are therefore unable to comment on how the ini¬ 
tial light curve extractio n process affects transiting plan¬ 
ets. We refer instead to IChristiansen et ^ (120131) for a 
discussion of pixel-level effects. They found that transits 
injected at the pixel level are usually recovered with high 
fidelity (final SNR = 96% — 98% expected SNR). 

The transit detection process as modeled in this paper 
consists of two distinct stages: (1) the putative event is 
identihed as a peak in the BLS periodogram and (2) the 
signal is accepted because a transit model provides a 5a 
improvement to a straight-line fit. We took advantage 
of the two-step nature of the search process when de¬ 
termining the search completeness for each star in our 
survey. 


For each star, we generated a set of 2000 trial plan¬ 
ets with orbital periods drawn from a log uniform dis¬ 
tribution extending from 0.5 to 200 days and uniformly 
distributed epochs of transit, radii (0.5 — 4 R 0 ), and im¬ 
pact param eters (0 — 1). We th en constructed tran¬ 
sit models (|Mandel fc Agoll I2002D for the trial plan¬ 
ets using the assigned planetary parameters and limb 
darkening parameters est imated from the coefhcients in 
iClaret fc Bloem^ (|201lD based on the stellar temper¬ 
atures and surface gravities. We resampled the transit 
models to the 29.4 minute long cadence integration time. 

The hrst question we wished to address was whether 
the transits of the trial planets would be accepted as 
5a detections in the ideal scenario in which the light 
curve was perfectly detrended and the orbital parame¬ 
ters were determined exactly. Accordingly, we first mul- 
tiplied the detrended light curves by the transit models 
(|Mandel fc Agoll[2002[) generated using the assigned trial 
planet parameters. We then checked whether the differ¬ 
ence in the between the best-fit transit model (which 
we guessed perfectly) and a no-transit model exceeded 
the 5cr detection threshold of = 30.863. If the tran¬ 
sit model was not preferred, then we recorded the trial 
planet as a non-detection. For a typical star, 8% of the 
trial planets were rejected at this stage. 

For the trial planets that would be accepted in the 
ideal scenario, we then conducted a more realistic test 
by multiplying the transit model by the raw PDC-MAP 
photometry before detrending usi ng the straight line ht 
method described in Section 13.21 We then re-checked 
whether the transit model was preferred at 5a. Trial 
planets that were not accepted at this stage were also 
recorded as non-detections. Overall, 4% of the trial plan¬ 
ets that were accepted in the ideal case of a perfectly 
detrended light curve were not accepted in this more re¬ 
alistic test. 

Finally, we tested whether the remaining trial plan- 
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ets would have been identified as peaks in the BLS pe- 
riodograms by running a full test for at least 25 trial 
planets for each star. We selected the trial planets for 
the full test by ranking the signals detected in the sec¬ 
ond round of testing in order of increasing where 

Ay^ was the value at which they were preferred to a non¬ 
transiting model. We chose the first 25 trial planets in 
the ranked list for which a random number draw yielded 
a result greater than 0.5. In other words, we thinned the 
sample of trial planets by 50% and selected those closest 
to the expected sensitivity threshold. If fewer than ten 
trial planets were recovered, we conducted up to 25 ad¬ 
ditional runs (for a total of 50 simulations) until at least 
ten planets were recovered. 

For the selected trial planets, we multiplied the corre¬ 
sponding transit model by the raw PDC-MAP photom¬ 
etry for the assigned host star and detrended the light 
curve using a 2000 minute median filter as explained in 
Section lXn We then fed the injected light curves into the 
detection pipeline described in Section [3] Although we 
considered the full range of planet periods (0.5-200 days), 
we halted the search process as soon as the injected signal 
was detected. If the signal was not detected, we termi¬ 
nated t he s earch using the usual conditions discussed in 

Section [321 



Fig. 5. — Empirical detection sensitivity versus the between 
fitting the detrended light curve with the injected transit model or 
with a no-transit model. The light blue histogram depicts the re¬ 
covery fraction for the 83699 BLS trial planets in bins of 1000 plan¬ 
ets. The solid red line marks the estimated likelihood of detec¬ 
tion for the 4398023 non-BLS injected planets predicted from the 
smoothed histogram of BLS results. 

6.1. Predicting Transit Detectability 

In total, we ran 83699 complete BLS injection simula¬ 
tions for the 2543 stars in our sample. We also injected 
604278 planets that had Ax^ below our 5a detection 
threshold. The remaining 4398023 injected planets had 
Ax^ above the detection threshold but were not tested in 
the full BLS simulation. We predicted the detectability 
of these trial planets by finding the fraction of BLS trial 
planets recovered as a function of the computed in 
the second round of transit model tests. We ranked the 
BLS trial planets by and computed the recovery 
fraction for each consecutive group of 500 planets. Next, 


TABLE 3 
Injected Planets 


KID 

Period 

(Days) 

Radius 

(R©) 

Insolation 

(F©) 

Recovery 

Status^ 

003835071 

0.686 

0.712 

409.813 

1.000 

005693298 

0.876 

0.531 

250.636 

1.000 

011560326 

1.900 

2.700 

48.871 

0.912 

008229458 

2.700 

0.788 

61.212 

0.520 

004931385 

20.420 

1.717 

3.056 

0.650 

009691776 

55.352 

2.555 

0.568 

0.904 

010032631 

85.665 

1.978 

0.664 

0.791 

002441562 

112.865 

0.925 

0.195 

0.000 

011013096 

155.207 

1.795 

0.123 

0.908 

002692704 

184.735 

1.181 

0.225 

0.000 


Note. — Table 0 is published in its entirety in the 
electronic edition of the Astrophysical Journal. A por¬ 
tion is shown here for guidance regarding its form and 
content. 

^ For the 83699 injected planets that were tested in the 
complete pipeline and the 604278 planets with be¬ 
low our 5cr detection threshold, the recovery status indi¬ 
cates whether the planet was detected (1 for recovered 
planets, 0 for unrecovered planets). For the remaining 
4398023 injected planets that had Ay^ above the de¬ 
tection threshold but were not tested in the full BLS 
simulation, the recovery status indicates the estimated 
likelihood of detection (see Section n. 

we smoothed the resulting histogram and predicted the 
likelihood of detection for the 4398023 non-BLS runs us¬ 
ing a cubic spline interpolation based on the smoothed 
histogram. We limited the maximum detection likeli¬ 
hood to 91.2%, which was the maximum value of the 
cubic spline in the histogram of the recovery rate for the 
full BLS runs. Figure |5| displays the histogram of the 
recovery rate for the BLS trial planets and the extrapo¬ 
lated likelihoods of detection for the non-BLS runs. 

6 .2. Assessing Pipeline Performance 

In total, we injected 5086000 transiting planets into the 
light curves of the 2543 stars in our sample. We provide 
a catalog of injected planet parameters and recovery re¬ 
sults in Table m Our pipeline successfully recovered 86% 
of signals injected with an expected SNR between 15 and 
20. For lower SNR, the pipeline performance decreased 
roughly linearly with anticipated SNR until reaching 52% 
recovery for signals with anticipated SNR between 5 and 
7. For the purpose of assessing pipeline performance as 
a function of SNR, we modeled the anticipated SNR of 
a transiting planet as: 

SNR = ——— Vntransit (2) 

KjiJr r transit 

where S is the median decrease in brightness during the 
injected transit, CDPPtransit is the Combined Differen¬ 
tial Photometric Precision (CDPP) on the timescale of 
a full transit of a planet on a circular orbit, and Utransit 
is the number of transits expected given the orbital pe¬ 
riod of the plan et and the number of day s the s tar was 
observed. As in iDressing fc Charbonneau! (|2013[) . we es¬ 
timated the CDPP on the timescale of a planetary transit 
by interpolating over the provided CDPP measured on 
3-, 6-, and 12-hour timescales. 

As a benefit of injecting multiple trial planets per star, 
we generated unique transit detectability maps for each 
star in our sample. For example, Figure [6] displays the 
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transit detectability maps for KID 7104554, a Kp = 15.3, 
Teff = 3957iP star with lower search completeness than 
the larger sample. We created the star-by-star transit 
detectability maps by gridding the injected planets in 
radius/period and radius/insolation space and calculat¬ 
ing the fraction of detectable planets within each grid 
cell. For the subset of trial planets that passed both 
tests yet were not selected for the full BLS search, 
we estimated the recovery fraction as explained in Sec¬ 
tion [O] 

After generating transit detectability maps in radius- 
period and radius-insolation space for each star indepen¬ 
dently, we created transit detectability maps for the full 
sample by summing the individual maps. Although the 
combined maps displayed in Figure [7] are useful for com¬ 
paring the sensitivity of any individual star to the sensi¬ 
tivity of the larger sample, the binning is rather coarse. 
We therefore generated a second set of combined sensi¬ 
tivity maps by sorting the full set of 5086000 injected 
planets into smaller grid cells in radius/period and ra¬ 
dius/insolation space. We then calculated the recovery 
fractions within each of the cells to produce the smoother 
sensitivity maps displayed in Figure [H 

As shown in Figure [51 we found that our pipeline is 
very sensitive to injected planets with radii larger than 
2.5 R 0 . Such planets were detected with nearly 90% ef¬ 
ficiency out to the maximum injected orbital period of 
200 days. Our pipeline had a significantly harder time 
detecting 1.0 — 1.5 R 0 planets with periods longer than 
100 days (recovery fraction = 36%) and 0.5—1.0 R 0 plan¬ 
ets with periods longer than 5 days (recovery fraction 
= 22%). Planets smaller than 1.0 R 0 were nearly un¬ 
detectable (recovery rate approximately 6 %) at orbital 
periods longer than 150 days. 

Inspecting the transit recovery map as a function of 
insolation revealed that the transit detectability changes 
sharply across the habitable zone (HZ). At the inner edge 
of the HZ (median orbital period of 50 days for the stars 
in our sample), we recovered 84% of 2.0R 0 planets and 
34% of 1.0 R 0 planets. At the outer edge of the habitable 
zone (median orbital period of 130 days), the sensitivity 
decreased to 80% for 2.0 R 0 planets and 25% for 1.0 R 0 
planets. This large change in sensitivity in a very inter¬ 
esting region of planet radius and insolation space reveals 
that the pipeline sensitivity within the habitable zone is 
not well described by a single number. 

6.3. Calculating Search Completeness 

The overall planet search completeness depends both 
on the detectability of a particular transiting planet and 
the likelihood that a particular planet will be observed 
to transit. We accounted for the latter factor by de¬ 
termining the mean geometric probability of transit for 
planets orbiting the stars in our sample at particular pe¬ 
riods or insolation levels. For a given orbital period, we 
computed the corresponding semimajor axis for a planet 
orbiting each of the stars in our sample. Next, we di¬ 
vided the stellar radii by the calculated semimajor axes 
to find the transit probability for a planet in a circular 
orbit. 

We then multiplied the transit probability by a correc¬ 
tion factor to account for the fact that the planets in our 
sample are more likely to b e eccentric and were accord - 
ingly more likely to transit (|Barnesll2007HKipDindl2014D . 


We adopted a cor rection fa c tor of 1.08 based on a beta 
distribution fit by iKinnind (I2013D to transiting planets 
with periods shorter than 382.3 days. Neglecting this cor¬ 
rection factor would lead to an underestimate of search 
completeness and an over estimate of th e planet occur¬ 
rence rate by roughly 8 % (jKiDT)indl2014[ l. The resulting 
search completeness plots are displayed in Figure [5] 

7. THE PLANET OCCURRENCE RATE 

In order to estimate the planet occurrence rate, we first 
generated smoothed maps of the detected planet popula¬ 
tion. For each planet candidate, we counted the number 
of links from the MCMC posteriors that fell within each 
grid cell in radius/period and radius/insolation spaceF^ 
When converting each link of the chains from light curve 
parameters to physical values, we accounted for uncer¬ 
tainties in the stellar parameters by drawing new stellar 
parameters from Gaussians centered at the reported val¬ 
ues with widths set by the reported errors. (In cases 
where the reported errors were asymmetric, we adopted 
the larger value.) We weighted each link so that the to¬ 
tal weight equaled one minus the false positive correction 
(see Section for a planet with the given radius. 

The errors on the planet radii and insolation flux were 
large enough that the posteriors from multiple candidates 
overlapped to produce smoothed distributions. For the 
orbital periods, however, the errors were small enough 
that each planet appeared isolated. For the purpose of 
calculating the planet occurrence rate, we artificially in¬ 
flated the spread of the period values so that the standard 
deviation of log^o R distribution was equal to 0.1. We 
then constructed smoothed distributions of the insolation 
flux received by each planet by converting the periods 
into semimajor axes using Kepler’s third law and stel¬ 
lar masses drawn from gaussian distributions centered 
on the reported value with widths set by the reported 
errors. 

As shown in Figure (TOl the detected planet popu¬ 
lation has peaks within the region P = 1 — 20 days 
and Rp = 0.7 — 2.5 R 0 . There is also a noticeable lack 
of large planets {Rp > 1.7R 0 ) and shorter periods 
(P < 2 days). In radius-insolation space (right panel 
of Figure [TUI), the highest peaks of the smoothed candi¬ 
date distribution are located at insolations of 20 — 50 and 
2.5 — 10 times the insolation received by the Earth. 

We estimated the planet occurrence rate by dividing 
the smoothed maps of the detected planet population 
in Figure [10] by the smoothed maps of the search com¬ 
pleteness in Figure [U] The resulting maps of the planet 
occurrence rate are displayed in Figure fTTI The division 
reveals that the lack of detected planets in the upper left 
corner of the left panel of Figure [TO] is quite meaning¬ 
ful. That region has very high search completeness, so 
the lack of detected planets in that region of parame¬ 
ter space implies that hot mini-Neptunes and Neptunes 
are rare around low-mass stars. At the opposite corner 
of the diagram in the small-planet, long-period regime, 
the relatively small number of detected planets does not 
indicate a low occurrence rate. On the contrary, those 

We di r ectly incorporated the posteriors for the KOIs fit by 
I Rowe et al.l ||2014I ~). For the other 15 KOIs with parameters drawn 
from previously published papers, we modeled the radii and periods 
by constructing Gaussian distributions using the reported values 
and errors. 
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Insolation 


Fig. 6.— Transit detectability maps for KID 7104554 as a function of planet radius and orbital period (Left) or insolation flux (Right) 
based on the results of our injection simulation. The small red plus symbols mark the 498 injected planets with below the 5c7 detection 
threshold. The 25 large circles indicate injected planets with above the detection threshold that were recovered (yellow, 17 planets) 
or undetected (red, 8 planets) during the full BLS test phase. The small gray plus symbols are the remaining 1477 injected planets with 
Ax^ above the detection threshold that were not selected for the full BLS test. The numbers within each cell denote the recovery fraction 
within the cell boundaries and the cells are color-coded so that darker colors correspond to lowe r detectability. The green dashed lines mark 
the maximum greenhouse (Max GH) and moist greenhouse (Moist GH) insolation limits fr omIKo pparapu et al.l ll2Q13b|j and the magenta 
dot-dashed lines mark the less conservative Recent Venus and Early Mars limits, also from fKoppm'am^^aLril2Q13bl') . 
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Fig. 7.— Combined transit detectability maps for the full stellar sample as a function of planet radius and orbital period (Left) or 
insolation flux (Right) based on the results of our injection simulation. The numbers within each cell denote the recovery fraction within 
the cell boundaries and the cells are color-coded so that darker colors correspond to lower detectability. These figures were produced by 
combining individual completeness maps for each star such as those displayed in Figure |6] As in Figure |6] the vertical lines in the right 
panel mark two definitions of the habitable zone. 


planets were detected despite relatively low search com¬ 
pleteness, so the underlying occurrence rate of such plan¬ 
ets is predicted to be high. 

Consulting the right panel of Figure [HI the estimated 
occurrence rate of small planets is highest at insolations 
below roughly 0.5 F 0 , but that is a region of low search 
completeness and the occurrence rate for such planets is 
not well constrained. In order to more readily see trends 
in the planet occurrence rate as a function of planetary 
properties, we binned the smoothed occurrence distribu¬ 
tions shown in FigurefTTlto produce the gridded diagrams 
shown in Figure I12I The gridded version of the radius- 
period diagram (left panel of Figure [T2ll clearly demon¬ 
strates that planet occurrence increases with decreasing 
planet radius and increasing logj^g P. 


7.1. Dependence on Planet Radius & Period 

Figure dsi Figure O and Table [4] display the planet 
occurrence rate as a function o f planet radius and orbital 
period. As in previous studi es (|Dressing fc Charbonneau! 
I2013t [Morton fc Swift]l201^ . we found that planets with 
radii > 3 R® are rare around small stars (at least out to 
orbital periods of 200 days). Concentrating on planets 
with periods shorter than 50 days, we observed a gen¬ 
eral trend of decreasing planet occurrence with increasing 
planet radius between 1 and 4 R^. There is an indi¬ 
cation in Figure [T4| that the planet occurrence rate may 
become flat in log^g P for periods longer than 10 days. 
However, the errors on the longest orbital bins are large 
enough that we cannot distinguish between a brief flat¬ 
tening between 10 — 100 days and a plateau extending 
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Fig. 8. — Smoothed maps of the fraction of injected planets that were detected by our pipeline. As indicated in the color bars, darker 
points correspond to lower detection fractions. Left: Planet radius versus period. Right: Planet radius versus insolation. As in Figure!^ 
the vertical lines mark two definitions of the habitable zone. 


Completeness Completeness 



-4.33 -3.42 -2.51 -1.61 -0.70 -3.69 -2.94 -2.18 -1.42 -0.67 


Fig. 9. — Smoothed maps of the search completeness accounting for both pipeline sensitivity and the geometric probability of transit. As 
indicated in the color bars, darker points correspond to lower search completeness. Left: Planet radius versus period. Right: Planet radius 
versus insolation. As in Figure the vertical lines mark two definitions of the habitable zone. 
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Fig. 10. — Smoothed distribution of planet candidates detected by our pipeline. The color scale is linear with lighter colors indicating a 
higher number of planets. Left: Planet radius versus period. Right: Planet radius versus insolation. The orange points with error bars 
are the planet candidates detected by our pipeline. As in Figure [gl the vertical lines mark two definitions of the habitable zone. 
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Fig. 11. — Smoothed plot of the derived planet occurrence rate as a function of planet radius versus orbital period (Left) or insolation 
(Right). Lighter colors indicate higher planet occurrence per grid cell and regions in which our pipeline detected < 15% of injected signals 
are marked in gray. As in Figure O the vertical lines in the right panel mark two definitions of the habitable zone. 
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Fig. 12. — Binned planet occurrence rate in period/planet radius space (Left) and insolation/planet radius space (Right). The numbers 
within each grid cell indicate the planet occurrence rate as a percentage (top) and the percentage of injected planets that were recovered 
by our pipeline (bottom). The gray regions have injected planet recovery rates below 15%. Some boxes have large Poisson errors; please 
see Tables HI-El As in Figure O the vertical lines in the right panel mark two definitions of the habitable zone. 
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Planet Radius (REarth) 


Fig. 13.— Planet occurrence (top) and cumulative planet 
occurrence (bottom) versus planet radius for planets with periods 
of 0.5 — 1.7 days (dark green), 1.7 — 5.5 days (teal), 5.5 — 18.2 days 
(light blue), 18.2 — 60.3 days (navy), and 60.3 — 200 days (pur¬ 
ple). The error bars are based on binomial statistics and the 
assume d smoo thing of the planet population. In this figure and 
Figures I14H15I we do not present occurrence rates for regions with 
pipeline sensitivity below 15%. 

out to much longer orbital periods. 

For orbital periods shorter than 50 days, we measure an 
occurrence rate of 0.56lo;o| Earth-size (1 — 1.5 R©) plan¬ 
ets and 0.461 o;q5 super-Earths (1.5 — 2R©) per small 
star. Extending the period range to 100 days, we es¬ 
timate 0.65to;o5 Earths and 0.57lQ;Qg super-Earths per 
star. Overall, we find 2.5 ± 0.2 planets per M dwarf 
with radii of 1 — 4 R© and periods shorter than 200 days. 
We provide cumulative planet occurrence rates for sev¬ 
eral additional choices of period and radius boundaries 
in Tabled 

7 . 2 . Dependence on Planet Radius & Insolation 

In Figure [15] and Tables |6l|7l we present the planet oc¬ 
currence rate as a function of stellar insolation and planet 
radius. In general, we find that planet occurrence in¬ 
creases both with decreasing planet radius (as discussed 
in Section EH) and with decreasing log^Q Fp. Intrigu- 
ingly, we observed that the occurrence rates of Earths, 
Super-Earths, and mini-Neptunes are comparable for in¬ 




Fig. 14.— Planet occurrence (top) and cumulative planet 
occurrence (bottom) versus orbital period for planets with radii 
ofO.5 —IRm (black), 1 —1.5Rm (dark gray), 1.5 —2.0Rm (brown), 
2 - 3 R® (o®range), and 3 - 4R® (red) 

solutions below roughly 30 F© but that planets larger 
than 1.5 R© were less common than smaller planets at 
insolations above 30 F©. The error bars on the coolest 
insolation bin are rather large, but the divergence of the 
Earth and Neptune occurrence relations might be due to 
photo-evaporation at short orbital periods. 

Across the size range we considered, the planet occur¬ 
rence rate versus log,Q Fp rises with decreasing insolation 
between 100—10 F© and appears roughly flat in log,Q Fp 
between 10 — 0.2 F©. Figure [TSl hints that the planet oc¬ 
currence rate might increase again at cooler insolations, 
but a larger sample of long period planets will be required 
to test that hypothesis. 

7 . 3 . The Occurrence of Potentially Habitable Planets 

The shaded regions in Figure [1^1 display one choice 
of habitable zone boundaries, but the definition of the 
“habitable zone” (HZ) is still rather uncertain. Tra¬ 
ditionally, astronomers have used the term to refer to 
the distance from the star at which li quid water could 
be present on the surfac e of a planet (|Dolel 11964 iHard 
119791 : lasting et H1I1993D . In theory, there could also be 
water-based life on worlds with subsurface oceans or non- 
water-based life on worlds like Titan, but any associated 
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TABLE 4 

Number of Planets Per Star Versus Orbital Period (In Percentage) 


(R-©) 

0.5 - 1.7 Days 

1.7 - 5.5 Days 

5.5 - 18.2 Days 

18.2 - 60.3 Days 

60.3 — 200 Days 

0.5 - 

1.0 

l-92lj-“l(70%) 

9-881^4^(54%) 

23.06llg9(33%) 

17.75l55gf (17%) 

(6%) 

1.0 - 

1.5 

1.83t°-®^(87%) 

9.1811^8(84%) 

< 26.54(79%)=* 

23.081158(64%) 

30.7019167(41%) 

1.5 - 

2.0 

0.33t^(90%) 

3.70llf^(88%) 

20.0611^5(87%) 

26.731199(84%) 

18.9011(1^(74%) 

2.0 - 

2.5 

< 0.13(90%) 

2.251157(90%) 

15.691191(89%) 

23.651151(88%) 

14.12liigf (84%) 

2.5 - 

3.0 

< 0.11(91%) 

1.1911^1(90%) 

6.541511(90%) 

10.42ll5|(89%) 

5.3011,55(88%) 

3.0 - 

3.5 

< 0.11(91%) 

0.58ll®l(91%) 

2.55ll«(91%) 

2.441111(90%) 

1.831152(89%) 

3.5 - 

4.0 

< 0.11(91%) 

0.36111^(91%) 

0.871155(91%) 

0.531151(90%) 

< 2.71(90%) 

0.5 - 

1.0 

1.38+";^^(66%) 

8.4211“ (44%) 

20.59ll^!^(26%) 

(11%) 

(3%) 

1.0 - 

1.5 

l-95tg-^?(86%) 

9.9411^3(83%) 

< 26.63(75%) 

26.851^1^9(58%) 

28.85l?l||(32%) 

1.5 - 

2.0 

0.4V°;®J(89%) 

4.15ll®|(88%) 

< 23.58(86%) 

24.5911^^(82%) 

19.98li54°7(67%) 

2.0 - 

2.5 

< 0.13(90%) 

2.7211^5(89%) 

18.7311^1(89%) 

27.5811^9(37%) 

18.081 J5^?i (82%) 

2.5 - 

3.0 

< 0.11(90%) 

1.591159(90%) 

8.2911,96(90%) 

14.5lll;g(89%) 

8.6111^^(87%) 

3.0 - 

3.5 

< 0.12(91%) 

0.65ll°°(91%) 

3.251137(90%) 

3-3711^(90%) 

1.9711,57(89%) 

3.5 - 

4.0 

< 0.11(91%) 

0.38im(91%) 

1.0511,59(91%) 

0.56ll?^(90%) 

< 2.34(89%) 

Note. — 

- In this table and all subsequent 

occurrence rate tables, the numbers in 

parentheses are the 


fraction of injected planets that were recovered within the given intervals. In addition, t he first set of entries 
are our estimates of the planet occurrence rate when using the stellar properties in the IHuber et al.l 1120141) 
catalog. The second set of entries (below the double line) are alternative e stimates constructed by revising 
the stellar radii to lie along an empirical temperature/radius relation from [Mann et all ll2Q13bl . see Section 
17.41 for details). 

^ We provide one-sigma upper limits instead of two-sided errors for grid cells with large Poisson errors and 
very low occurrence rates. 


TABLE 5 

Cumulative Number of Planets Per Star Versus Orbital Period (In Percentage) 


7?p(R©) 0.5 — 10 Days 0.5 — 50 Days 0.5 — 100 Days 0.5 — 150 Days 0.5 — 200 Days 


0.5 

- 1.0 

22.08111^(54%) 

51.4311,99(44%) 

54.73111^(40%) 

59.71111^(37%) 

59.711144(36%) 

1.0 

- 1.5 

21.18llfj(84%) 

56.251151(79%) 

65.091152(75%) 

82.17ll°^(72%) 

88.74ll,5|(70%) 

1.5 

- 2.0 

10.76lljg(89%) 

45.95111^^(87%) 

57.3511^^(86%) 

66.581111(85%) 

69.72im(84%) 

2.0 

- 2.5 

8.681191(90%) 

35.62ll°9(89%) 

49.811115(89%) 

54.591146(88%) 

55.721152(88%) 

2.5 

- 3.0 

3.89l9;|g(90%) 

15.4611^^(90%) 

21.9711,57(90%) 

23.19llJ4(90%) 

23.451199(90%) 

3.0 

- 3.5 

1.731154(91%) 

4.651119(90%) 

7.151159(90%) 

7.381144(90%) 

7-4111,6? (90%) 

3.5 

-4.0 

0.721198(91%) 

1-6311^1(91%) 

2.1llgj^(91%) 

2.2711^1(90%) 

2-2711,0? (90%) 


0.5 

- 1.0 

18.3911?^ (49%) 

39.9311?? (38%) 

40.601?;?? (34%) 

40.6011?? (31%) 

40.6011?? (30%) 

1.0 

- 1.5 

22.8511?? (82%) 

60.7611?? (75%) 

68.38ll?|(71%) 

85.73l?;?|(68%) 

91.47 ± 5.04(66%) 

1.5 

- 2.0 

11.96iy?(88%) 

45.8711?? (86%) 

57.261165(85%) 

66.74l?;®5(83%) 

70.2311?? (82%) 

2.0 

- 2.5 

9.98114? (90%) 

42.481?;?? (89%) 

58.49l?-9?(88%) 

65.45ll5?(88%) 

67.1211?? (87%) 

2.5 

- 3.0 

4.94l?;5|(90%) 

20.6015;?? (90%) 

29.8511?? (90%) 

32.4011?? (89%) 

33.0111?? (89%) 

3.0 

- 3.5 

2-0911,?? (91%) 

6.321?;?? (90%) 

8.671?;?? (90%) 

9.131?;?? (90%) 

9.251?;?? (90%) 

3.5 

-4.0 

0.83lli?(91%) 

1.92ll?5(9i%) 

2.1411,?? (90%) 

2.20l?;??(90%) 

2.2ll?;??(90%) 


Note. — As in Table [4l th e entries below the double horizontal line are the estimates 1 
stellar radii (see Section l7.4ll . 


biosignatures would be difficult to interpret remotely. 
Accordingly, astronomers have concentrated thus far on 
the search for life as we know it, meaning surface-based 
lifeforms that depend on liquid water and might gen¬ 
erate biosignatures that could alter the composition of 
their homeworld’s atmosphere. 

Even within that rather narrow definition, there are 
many assumptions that can affect the choice of habit¬ 
able zone boundaries. In particular, the assumed mass 


and composition of the planet’s atmosphere affects the 
surface pressure and therefore the t emperature range a t 
which water would be liquid (e.g.. IVladilo et al.ll2013ll . 
For instance, iPierrehumbert fc Gaidos| (120111) showed 
that planets with thick hydrogen atmospheres would 
have sufficient surface pressure to retain surface liquid 
water out to distances of 2.4 AU. It is uncertain whether 
biosignatures would be detectable i n such atmospheres 
(|Seager et al.l 120131 iHu et al.l I2013I1 , but those worlds 
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TABLE 6 

Number of Planets Per Star Versus Insolation (In Percentage) 


(n©) 

0.2 - 1.1 F® 

1.1 -6.3F® 

6.3 - 35.6 F® 

35.6 - 200 F® 

0.5 

- 1.0 

(11%) 

22.63l“6f(24%) 

23.481®-^® (44%) 

3.52ll®J(60%) 

1.0 

- 1.5 

24.02(47%) 

< 35.71(69%) 

13.45l2-®g(81%) 

5.36ll®®(85%) 

1.5 

- 2.0 

17.70tJ\™(75%) 

< 31.17(85%) 

10.5913-^1(88%) 

1.38ll®3(89%) 

2.0 

- 2.5 

13.97+®'®^ (85%) 

20.16t®®°(88%) 

9-8413-34(89%) 

0.751131(90%) 

2.5 

- 3.0 

5.47t®;®^(88%) 

8.90l®;°^(90%) 

3.941128(90%) 

0.42ll®2(90%) 

3.0 

- 3.5 

1.7lt^'^®(89%) 

2.18t®-J®(90%) 

1.7llJ-®®(91%) 

0.27llf3(9i%) 

3.5 

-4.0 

< 2.33(90%) 

0.8012-28(91%) 

0.5711-11(91%) 

0.19l°-ii(91%) 

0.2 

- 1.1 

1.1-6.3 

6.3-35.6 

35.6-200.0 


0.5 

- 1.0 

(6%) 

10.46l®®®(16%) 

17.22l®-i°(32%) 

3.23ll^®(52%) 

1.0 

- 1.5 

26.87j:g®g®’’(37%) 

29.27l®12(g3%) 

13.19ll^®(78%) 

5.8411^® (84%) 

1.5 

- 2.0 

18.99tJ®g°®(71%) 

24.541^-3^(83%) 

11.35ll®3(87%) 

l-79ll®I(88%) 

2.0 

- 2.5 

15.29(ti°4" (83%) 

22.531^-21(88%) 

10.9311®® (89%) 

1.05l°-®i(90%) 

2.5 

- 3.0 

7.68^3 36(87%) 

11.831®-^® (89%) 

4.40111® (90%) 

0.4711®® (90%) 

3.0 

- 3.5 

l-90to8?(89%) 

2-7013-3^(90%) 

2.1111-33(90%) 

0.20+11^(91%) 

3.5 

-4.0 

< 1.81(90%) 

0.76l2;22(90%) 

0-75112^(91%) 

0-1911,1® (91%) 


Note. — As in Table |4] the entries below the d ouble horizontal line are the estimates 
based on the revised stellar radii (see Section |7.4| |. 


TABLE 7 

Cumulative Number of Planets Per Star Versus Insolation (In Percentage) 


itp ( R® ) 

0.2 - 200 F® 

1.0- 200 F® 

10.0 - 200 F® 

50.0 - 200 F® 

100.0 - 200 F® 

150.0 - 200 F® 

0.5 

- 1.0 

63.36llj^(34%) 

49.6311,“ (41%) 

17.7911,®^ (51%) 

2.09111^(58%) 

0.6i11«!(61%) 

0.061°,?® (62%) 

1.0 

- 1.5 

74.671111(70%) 

50.651®;!® (77%) 

15.00111^(83%) 

3.5011,33(85%) 

1.00ll«^(85%) 

0.16l°-,3?(86%) 

1.5 

- 2.0 

57.6lll®^(84%) 

39.921®;^® (87%) 

6.60lljJ(88%) 

0.901°,^! (89%) 

0.2311?? (89%) 

0.0413-23(89%) 

2.0 

- 2.5 

44.7211®® (88%) 

30.74ll“(89%) 

6.41111^(90%) 

0.33111^(90%) 

0.0211 J® (90%) 

0.001°-,°® (90%) 

2.5 

- 3.0 

18.72ll®2(90%) 

13.26ll®2(90%) 

2.85111! (90%) 

0.2311!® (90%) 

O.Olll®? (90%) 

< 0.00(90%) 

3.0 

- 3.5 

5.8711^^(90%) 

4.1611^2 (9o%( 

1.2811,23(91%) 

0.19ll!®(91%) 

0.0211 J® (91%) 

< 0.00(91%) 

3.5 

-4.0 

2.031121(90%) 

1.5611^3(9^%) 

0.49ll^®(91%) 

0.18ll!3(91%) 

< 0.00(91%) 

< 0.00(91%) 

0.5 

- 1.0 

34.0811,31(28%) 

30.911111^(34%) 

15.1311,'® (44%) 

2.03llli^(52%) 

0.5411®! (55%) 

0.051?;?? (57%) 

1.0 

- 1.5 

75.17ll,®J(64%) 

48.30llg2(73%) 

15.6111,^^(81%) 

3.95113® (83%) 

1.30111? (84%) 

0.22l?-??(84%) 

1.5 

- 2.0 

56.6611^^(81%) 

37.67llgg(85%) 

7.261111(88%) 

1.141134(88%) 

0.2811?? (88%) 

0.041°;?? (88%) 

2.0 

- 2.5 

49.80ll^®(87%) 

34.5111^1(89%) 

7.3i11®3(89%) 

0.5411®® (90%) 

0.0211“ (90%) 

< 0.00(90%) 

2.5 

- 3.0 

24.38llj^(89%) 

16.7111^^(90%) 

3.17ll,?|(90%) 

0.35111^(90%) 

0.0511?? (90%) 

< 0.00(90%) 

3.0 

- 3.5 

6-9111.36(90%) 

5.0111^1(90%) 

1.5i11,^?(91%) 

0.17ll,!3(91%) 

0.0311?? (91%) 

< 0.00(91%) 

3.5 

-4.0 

1.9311,33(90%) 

1.69ll3|(91%) 

0.62l°-®2(91%) 

0.18ll,!!(91%) 

< 0.01(91%) 

0.001°;°® (91%) 


Note. — As in TableU the entries below the double horizontal line are the estimates based on the revised stellar radii 
(see Section [7a . 


could still be habitable. 

The presence of clouds adds an additional complication 
by both cooling and heating the planet. Clouds are par¬ 
ticularly important in the case of tidally-locked planets, 
which might be a common fate for pla nets orbiting within 
the habitable zones of M dwarfs. I Yang et al.l 1)20130 
demonstrated that a tidally-locked M dwarf planet might 
develop a persistent cloud patch above the sub-stellar 
point. That cloud patch would have a higher albedo 
than the planetary surface and would allow the planet to 
be much cooler at a given separation than a cloud-free 
model would predict. As a result, the habitable zone 
for a cloudy, tidally-locked planet could extend to inso¬ 


lations as high as Fp = 1.76F0 for the moist greenhouse 
limit rather than t he limi t of Fp < 0.88 F® calculated by 
iKopparapu et al.l (|2013al) for a cloud-free model. Even 
for non-tidally-locked planets, the presence of clouds can 
expand the distances corresponding to the boundaries 
of the habitable zone by roughly 40% d epending on the 
degree of cloud coyer (iSelsis et al.ll2007[) . 

The orbital geometry of exoplanets is also impor¬ 
tant when assessing planetary habitability. For in¬ 
stance, planets with high obliquities or eccentric or¬ 
bits might be partially habitable at certain latitudes 
or during certain times of y ear (iWUliams^^KAstingi 

119971 : IWilliams fc Pollardll200^ iSoiegel et al.ll2008ll2009l : 
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Dressing & Charbonneau 


Dressiiig_et_ab I l201Ct iCowan et al.l 120121: iDobroyolskisI 


20131 : lArmstrong et al.l 12014 Linsenmeier et alJ 120151) . 
Depending on the timescale for the temperature of the 
planet to change (which depends on factors such as the 
fraction of surface covered by ocean), such planets may 
undergo periodic global glaciations punctuated by short¬ 
lived epochs d uring which the surfa c e is warm enough fo r 
liquid water (iPierrehumbertI [20051 : iSoiegel et al.ll201^ . 
In addition, the primordial obliquities of close-in planets 
orbiting M dwarfs may be significantly eroded by tides 
(jHeller et al.ll20lTh . 

Constructing a multi-dimensional habitable zone 
model for planets orbiting M dwarfs is beyond the scope 
of this paper, but we aspire to provide enough informa¬ 
tion so that other researchers can assess the abundance 
of planets within their chosen habitable zone boundaries. 
We therefore provide occurrence rates for a few possible 
choices of habitable zone boundaries in Table [51 

In the most conservative case, we adopt the maximum 
greenhouse (Max GH) and moist greenho u se (Mo ist GH) 
insolation limits from iKopparapu et all (|2013b[l . The 
Max GH limit is the insolation at which adding addi¬ 
tional CO 2 can no longer heat the surface of the planet 
because Rayleigh scattering begins to dominate over the 
greenhouse effect. At the inner edge, the Moist GH 
limit corresponds to the insolation at which the planet’s 
stratosphere becomes dominated by water vapor. At that 
point, the planet’s reservoir of hydrogen quickly escapes 
to space. 

Table [8] also provides estimates based on the assump¬ 
tion that Venus and Mars were habitable at earlier times 
in their histories. For a Sun-like star, those constraints 
correspond t o insolation limits of 1.77 6 Fm and 0.321 Fq, 
respectively (jKopparapu et al.ll2013bllal ). The insolation 
boundaries are lower for planets orbiting M dwarfs be¬ 
cause the incoming radiation is redder. For a typical 
star in our sample {T^s = 3748K), the boundaries are 
1.543 F 0 and 0.228 Fq, respectively. 

Even more optimistically. Table [8] include s HZ occur¬ 
rence rates using the cloudy inner HZ from lYang et al.l 
dMl. tw o cho ices of desert world albedos from 
IZsom et al.l (I2013fl. and the conven ient limits of 0.25 — 
4F0 used bv iPetigura et al.l (|2013all . We do not provide 
an estimate based on the hyd rogen atmosphere HZ of 
iPierrehumbert fc Gaidod (|2011ll because our search com¬ 
pleteness is very low at the maximum allowed separation 
of 2.4 AU. 

The appropriate radius range to consider for a po¬ 
tentially habitable planet is more clearly defined than 
the appropriate insolation range. Based on radial veloc- 
ity follow-up observations of Kepler planet candidates, 
IRogersI (120151 ) argued that the majority of planets larger 
than 1.6 R 0 contain too many volatiles to be rocky. 
This result agrees with pre vious fits to me a sured exo¬ 
planet masses and radii bv iWeiss fc Mar"^ (|2014ll and 
simulations by iLopez fc: FortneA^ 20141) ^ E^urthermore, 
[Dressing et al.l ( 201511 found that all five exoplanets 
smaller than 1.6 R 0 with masses and radii measured to a 
precision better than 20% have densities consiste nt with 
an Ea rth-like mixture of iron and silicates. Like [Rogers! 
(|2015f) . they noted that planets larger than 1.6R 0 radii 
have densities inconsistent with rocky compositions. Due 
to observational constraints, the population of planets 
with well-constrained densities is strongly biased towards 


highly irradiated planets. We therefore include a broader 
range of radius choices in Table [ 8 ] to account for the pos¬ 
sibility that the transition between rocky and gaseous 
planets might occur at a slightly different radius for less 
irradiated planets. For instance, the 2.35 R 0 exoplanet 
Kepler-10c has a measured mass of 17.2 ± 1.9 M 0 and a 
bulk density of 7.1 ± lgc m~^, higher than the d ensities 
of most 2 — 3 R 0 planets (iDumusaue et al.l[2014[l . 

Adopting the most conservative assumptions 
(1.0R 0 < Rp < 1.5 R 0 , outer HZ = Max GH, in¬ 
ner HZ = Moist GH), we estimate an occurrence rate 
ofo.i 6 l°:J^ potentially habitable 1 — 1.5 R 0 planets per 
M dwarf. The predicted occurrence rates of super-Earths 
(1.5 — 2.0R 0 ) and larger planets ( 2 R 0 < Rp < 4 R 0 ) 
within the same habitable zone boundaries are 0 . 12 ^q q 5 
super-Earths and O.lb^ons larger planets per M dwarf. 
Expanding the radius range to 1 — 2 R 0 or increasing 
the habitable zone boundaries to the limits for recent 
Venus and early Mars increases the assumed occurrence 
rate to 0.27to'o8 O.24lg0|, respectively. In the 

most optimistic case, we estimate an occurrence rate of 
0 . 86 lg'Q 4 desert worlds with albedos of 0.8 and radii of 
1 — 2 Rffi receiving insolations between the IZsom et al.l 
(|2013f) inner limit and the Max GH outer limit. The as¬ 
sumed occurrence rate of potentially habitable M dwarf 
planets therefore varies by a factor of four depending on 
the specific choice of radius and insolation boundaries. 

The range of HZ possibilities in Table |5| is particu¬ 
larly useful for compar ing our results to those of previous 
studies. For instance, iPetigura et al.l (|2013all estimated 
that 22% of FGK stars host 1 — 2 R 0 planets receiving 
0.25 — 4 F 0 . Within the same boundaries, we find an 
occurrence rate of 0.83to o 5 small planets per M dwarf 
HZ. 

The difference in our estimates might suggest that hab¬ 
itable zone pl anets are more c ommon around lower-mass 
stars, but the IPetigura et all ()2013all prediction is based 
on an extrapolation of the occurrence rate for shorter pe¬ 
riod planets to longer period orbits assuming that planet 
occurrence is flat in logP. If the planet occurrence rate 
actually increases with log P at longer periods, then per¬ 
haps the occurrence rates of potentially habitable plan¬ 
ets orbiting FGK and M dwarfs are more similar. Such 
a change in the slope of the FGK star planet occur¬ 
rence rate at longer periods could be explained by the 
radial- and temperature-dependence of the physics gov¬ 
erning planet formation. 


7.4. Implications of Systematic Biases in Modeled 
Stellar Radii 


The stellar parameters for the majority of the stars 
in our sample were estimat ed by fitting Dartmouth 
stellar models to photometric ((Dressing fc Gharbonnead 
(2011: [Gaido^[2nil iHuber et al.ll201^ or spectroscopic 
( Mann et al.l 20lS ^ lMuirhead et al.ll2012al) observations. 
The except i ons ar e one star with parameters from 
iMann et al.l (l2013bll and two very-lo w mas s stars with 
parame ters f rom Martin et al.l (I2013D . Both IMann et al.l 
(|2013bD and iMartfn et al.l ~ 2013!) estimated stellar radii 


using empirical relations based on interferomet ric obser¬ 
vations of low-mass stars (iBovaiia n et al.ll 2012f). 

Several recent s t udies fe.g.. iBovaiian et al.l 120121 : 
IMann et al.l l2013aL iNewton et al.l 1201511 have demon- 
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TABLE 8 

Habitable Zone Occurrence Rates (In Percentage) 


Fp{F(s) 

0.25 - 0.88 

0.23 - 1.54 

0.25-4.00 

0.25 - 1.76 

0.25 - 2.78 

0.25 - 5.85 

Outer HZ: 

Max GH*" 

Early Mars'^ 

Fixed® 

Max GH*’- 

Max GH^^ 

Max GH*’- 

Inner HZ: 

Moist GH'" 

Recent Venus^ 

Fixed® 

Gloudy Moist GH‘= 

Desert (a=0.2)4 

Desert (a=0.8)‘l 

0.5 - 

1.0 R0 

(13%) 

(14%) 

25.28/1/00(18%) 

15.75/1*4^(15%) 

20 . 61 / 1 * 3 !® (17%) 

27.56/1®87(20%) 

0.8 - 

1.0 R® 


13.09/1^7^(22%) 

20.95/1®4?1(28%) 

13.30/1*77(24%) 

17-23/1/11(26%) 

20.05/1®97(30%) 

1.0- 

1.5 R® 

15 . 82 / 1 ^ 5 “ (48%) 

24.28/1/®® (50%) 

46.77/1®42®(56%) 

21.82/1*57(52%) 

31.65/7o7(55%) 

46.90/1/7(59%) 

1.5 - 

2.0 R® 

11.54/®;^^(75%) 

20.69/g°3®°(76%) 

36.07/1°9°0(79%) 

21.48/1°27(77%) 

28.04/10g7(78%) 

39.35/1;7(80%) 

2.0- 

2.5 R® 

10 . 25 /®;®® (85%) 

17.09///® (85%) 

29.17//33(86%) 

18.23//11(86%) 

23.99///1(86%) 

30.97/1;7(87%) 

2.5 - 

4.0 R® 

5 . 03 / 1 ^? (89%) 

10.58//J® (89%) 

16.6l//^l(90%) 

11.09//?1(89%) 

13.79/1;11(89%) 

17.96/1;|1(90%) 

1.0- 

2.0 R® 

27.36/1^4^(61%) 

44 . 97 / 9 * 2 ^( 63 %) 

82.84//®® (68%) 

43.31/1/7(64%) 

59.68/l®9l®(67%) 

86.25/l;ll(69%) 

2.0- 

3.0 R® 

13.88/®;^^ (86%) 

24 . 44 / 105 ^( 87 %) 

41.17//®1(87%) 

25.97/l°5li(87%) 

33.87/10g7(87%) 

44.16/1;11(88%) 

3.0- 

4.0 R® 

1.40//^® (90%) 

3 . 23 //°® (90%) 

4.6l//?l(90%) 

3.34//7(90%) 

3.91//7(90%) 

4.77//7(90%) 

2.0 - 

4.0 R® 

15.28/®;g(88%) 

27.67/J°8^'‘(88%) 

45.78/l°o°®(89%) 

29.32/l°g®i(88%) 

37.78/797(89%) 

48.93/l;|3(89%) 

0.5 - 

1.4 R® 

19.85/i®7^%28%) 

37 . 34 /®/®® (30%) 

67.77/1®9?®(35%) 

35.60/l®g7(31%) 

49.48/1*37(33%) 

69.67/1/7(37%) 

0.5 - 

1.0 R® 

(8%) 

(9%) 

(12%) 

(10%) 

(11%) 

(14%) 

0.8- 

1.0 R® 

(14%) 

4 . 49 / 109 ®! (15%) 

8 . 55 / 1191 '! (20%) 

5.54/757(16%) 

7.76/1/7(18%) 

9.12/1/7(22%) 

1.0 - 

1.5 R® 

19.97/?°9f (38%) 

28.20/®igf (41%) 

49.46/1/10(48%) 

26.25/1®87(43%) 

35.71/1*91® (46%) 

48.68/1®37(51%) 

1.5 - 

2.0 R® 

13.17/1/^(69%) 

20 . 28 /^ 15 ^( 71 %) 

34.10/l°4f (75%) 

20.32/737(72%) 

26.12/1070(74%) 

36.79/1;7(76%) 

2.0- 

2.5 R® 

11.62//®® (83%) 

18.38/^o^g®(83%) 

32.50//?^(85%) 

19.41/®;®! (84%) 

26.23/0-01(84%) 

34.48/®-/! (85%) 

2.5 - 

4.0 R® 

7.07///® (88%) 

12.65///^ (89%) 

21.53//^! (89%) 

13.42//®! (89%) 

17.56/®;®® (89%) 

23.29/^;7(89%) 

1.0 - 

2.0 R® 

33.14/i®6f(54%) 

48.48/1/®® (56%) 

83.55//l®(61%) 

46.57/757(57%) 

61.83/1®47(60%) 

85.47//jl(63%) 

2.0 - 

3.0 R® 

17.14/i°8f (85%) 

28.28/liif (85%) 

49.13/l°2f (86%) 

29.93/1/7(86%) 

39.87/737(86%) 

52.58/®;!® (87%) 

3.0 - 

4.0 R® 

1.54//^2(89%) 

2.75/®'/! (89%) 

4.89//08(9o%) 

2.89/®;7(89%) 

3.92/®;079o%) 

5.18//o®(90%) 

2.0- 

4.0 R® 

18.69/g/f (87%) 

31.03/l.i40'i(87%) 

54.03/l°2l®(88%) 

32.83/737(88%) 

43.79/747(88%) 

57.77/1;7(88%) 

0.5 - 

1.4 R® 

18.48/®°5f (21%) 

29 . 91 / 9 / 10 ( 22 %) 

53.10/1/7(27%) 

28.96/1®37(24%) 

39.97/1*57(26%) 

53.79/!®57(30%) 


Note. — As in Table[^ the entries below the double horizontal line are the estimates based on the revised stellar radii (see Section l7.4ll . 

^ These habitable zone limits are from fKopparapu et all ll2013bD . 

^ These are the “simple” habitable zone boundaries adopted by Petigura et al.l ll2013ah . 

^ This limit approximates the effect of clouds iYang et al. 12014 ) by increasing the flux at the inner edge of the habitable zone a factor of two 
compared to the baseline calculation bv IKopparapu et al.l i 2013bD . 

^ For these calculations, the inner edge of the habitable zone was set to the values predicted bv IZsom et al ] l[20ll) for hot “desert” worlds with 
low relative humidity. 


strated that theoretical stellar models do not accurately 
reproduce the observed radii o f low-mass stars. As ex¬ 
plained in iNewton et al.l (|2015[ 1 , there are two main is¬ 
sues: 

1. The radii of model stars with T^ff < 4000K are 
smaller than the interferometrically-measured radii 
by approximately 0.04 — 0.09 Rq. 

2. Variations in metallicity produce significant 
changes in the modeled radii of low-mass stars 
whereas observations reveal that metallicity actu¬ 
ally has little influence on the radii of low-mass 
stars. 

Although the Dartmouth stellar models perform better 
than many alternative models, both of these effects may 
have caused the radii of the stars in our sample to be 
systematically underestimated. In order to gauge the 
magnitude of this effect, we recalculated the radii for our 
stellar sample using an empirical temperature/radius re¬ 
lation for main sequence stars with 3300K < T^ff (Equa¬ 


tion 6 in IMann et al.l I2013al with the additio nal signif¬ 
icant figures reported bv INewton et al.l 120151) . We did 
not consider c hange s in the stellar temperatures, but 
INewton et al.l (|2015f) demonstrated that the tempera - 
tures we estimated in iDressing fc Charbonneaul (|2013[ 1 
were consistent with predictions based on empirical ob¬ 
servations (our values were lower by 40 ± IlOK). Using 
the empirical temperature/radius relation to revise the 
radii of the 2437 stars in our sample with T^s > 3300, 
we found that the median change in radius (Ai?*) was 
an increase of 0.026 R© (6%). The change was highly de¬ 
pendent on the assumed metallicity; stars with assigned 
[Fe/H]< —0.5 displayed a median size increase of 0.05 R© 
(ll%) while the estimated radii of stars with assigned 
[Fe/H] > 0 shrank by 0.016R© (5%). 

For the planet host stars in our sample, the in¬ 
crease in the stellar radii leads to larger predicted radii 
and increased insolation fluxes for the associated planet 
candidates. The median planet radius increase was 
6 .6%, but the amplitude of the change varied consid- 
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Fig. 15. — Planet occurrence (top) and cumulative planet occur¬ 
rence (bottom) versus insolation for planets with radii of 0.5 — 1 R 0 
(crimson), 1 — 1. 5 R 0 (orange), 1.5—2.0Rm (brown), 2 — 3 R 0 (dark 
gray), and 3 —4R0 (black). As in Figure!^ the vertical lines mark 
two definitions of the habitable zone. 

erably. The systems most strongly affected by the re- 
vision of the stellar radii were: KOI 3102 (+48%), 
KOI 2650 (+26%), KOI 2418 (+20%), KOI 2006 (+17%), 
and KOI ^2 (+16%). Two of the KOIs in these sys¬ 
tems (KOI 3102.01 and KOI 2650.02) were missed by 
our planet detection pipeline so they did not enter into 
our calculation of the planet occurrence rate. 

In addition to altering the radius estimates for the de¬ 
tected planet candidates, the changes in the stellar radii 
affect the estimated survey completeness and, in turn, 
the derived occurrence rate. If the stellar radii are typ¬ 
ically 6% larger, then the search completeness we dis¬ 
played in Figure [7] for planets between 0.50 — 4.00 R 0 
actually corresponds to 0.53 — 4.24 R® planets. We ac¬ 
counted for this effect by generating new search com¬ 
pleteness maps following the procedure outlined in Sec¬ 
tion [6] but correcting the radii and insolation flux envi¬ 
ronments of the injected planets to reflect the new radius 
estimates for each star. We then recalculated the planet 
occurrence rate using the updated search completeness 
maps and the revised planet properties. We present the 
corresponding planet occurrence maps in Figure [16] and 
include the resulting planet occurrence rates below the 


double horizontal lines in Tables [HIHl 

As expected, the most noticeable difference between 
the occurrence maps displayed in Figures [T2| and [T6| is 
that the ridge of high planet occurrence has moved up¬ 
ward to large radii. Similarly, the region of low search 
completeness now encompasses a slightly larger portion 
of our chosen parameter space. Using the revised stellar 
radii, we calculated occurrence rates of 0.611 qq 5 Earth- 
size planets and 0.461 q q 5 super-Earths per low-mass 
star with periods shorter than 50 days. These rates are 
nearly identical to the estimates presented in Section 0 
Within the habitable zone, we estimate a frequency of 
0.201003 Earths and 0.13lg;j!)5 super-Earths per star 
when adopting the moist G H inner limit and t he max - 
imum GH outer limit from iKopnarapu et al.l (l20I3bll . 
These estimates are 26% and 14% higher, respectively, 
than the rates of O.ldlgy Earths and O.12lg05 super- 
Earths per HZ presented in Section |7l Although increas¬ 
ing the assumed stellar radii alters the inferred occur¬ 
rence rates, the dominant source of error is the relatively 
small number of potentially habitable small planets. 

8. SUMMARY & CONCLUSIONS 

In this paper, we presented an updated estimate of 
the planet occurrence rate for early M dwarfs based on 
the full four-year Kepler ddXa. set. We developed our 
own planet detection pipeline to search for transiting 
planets in the Aep/er light curves. We then character¬ 
ized the completeness of our pipeline by injecting simu¬ 
lated transiting planets into the Kepler light curves and 
attempting to recover them. Our search of the light 
curves of 2543 small stars with at least 1000 days of Ke¬ 
pler photoTaetry revealed 3215 possible planetary tran¬ 
sits. We thoroughly inspected all available follow-up 
observations for these objects and accounted for transit 
depth dilution for systems with close stellar companions. 
We accepted 156 planet candidates, one of which was not 
a previously known Kepler planet candidate. 

We then measured the occurrence rate of small plan¬ 
ets around small stars by dividing smoothed maps of 
the detected planet population by maps of our pipeline 
search completeness in radius-period space and radius- 
insolation space. We found that Earth-sized planets 
(1.0— 1.5 R®) are common and calculated an occurrence 
rate of 0.56tg g 5 Earth-sized planets with periods shorter 
than 50 days per early M dwarf. We also found an oc¬ 
currence of 0.46^005 super-Earths (1.5 — 2 R 0 ) with pe¬ 
riods shorter than 50 days per early M dwarf. For or¬ 
bital periods shorter than 200 days and planet radii of 
1 — 4 R 0 , we estimated a cumulative planet occurrence 
rate of 2.5 ± 0.2 planets per M dwarf. 

Within a conservatively defined habitable zone based 
on the moist greenhouse and maximum greenhouse lim¬ 
its ([KoDDaraDU et al.l 1201310 we estimated occurrence 
rates of O.16lg07 Earth-size (1.0 — 1.5 R 0 ) planets and 
O.12lg05 (1.5 — 2.0R 0 ) super-Earths per small star. 
Adopting a wider planet size range of 1 — 2 Rm and con¬ 
sidering the effects of clouds (lYang et al.ll2013fl increased 
our estimate to 0.43^Q gg potentially habitable planets 
per star. Gonsidering desert worlds (|Zsom et al.ll2013ll 
would increase the measured occurrence rate to nearly 
one potentially habitable planet per M dwarf. These 
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Fig. 16. — Alternative calculation of the planet occu rrence rate in period/planet radius space ( Left) and insolatio^planet radius space 
(Right) using the revised stellar radii (see Section l7.4|l . The annotations are the same as in FigureAs in Figurethe vertical lines in 
the right panel mark two definitions of the habitable zone. 

estimates span the range of previous estimates of the oc¬ 
currence rate of potentially habitable M dwarf planets. 

An order of magnitude calculation multiplying the oc¬ 
currence rate of potentially habitable 1 — 1.5 R® planets 
between the empirical early Mars outer boundary and 
the recent Venus inner boundary by an estimate of the 
number density of small stars in the galaxy from the 
RECONS survey (|Henrv et al.ll2006tlWinters et al.ll2015ll 
therefore suggests that the nearest potentially habitable 
planet is most likely 2.6 ± 0.4 pc away and is within 
3.5 pc with 95% confidence. This estimate assumes that 
the occurrence rate of potentially habitable planets orbit¬ 
ing later M dwarfs is identical to that for early M dwarfs , 
which is consistent with the results of iBerta et al.l (|2013D . 

Correcting for the geometric probability of transit 


(assuming that 1.4% of potentially habitable M dwarf 
planets transit), the nearest transiting potentially 
habitable planet is likely to be 10.6l( g pc away and is 
within 14.6 pc with 95% confidence FI Early M dwarfs 
at distances of 3.5 pc and 14.6 pc would have apparent 
K band magnitudes of 2.9 and 6.0, respectively, well 
within the magnitude range probed by current and 
upcoming planet surveys of nearby, brigh t stars such 
as CARMENES jO^enbach et all I^HOl. CHEOP S 
(jBroeg et al.| |201 3[i. Exopla netSat (iSmith e t al.l 120 1011 
ExTr A (jBonfils et al.l [2f)l^. HPF (iMahadeyan et al.l 
201011 ■ MEarth (Nutzman & Charbonneau 120081 : 
Berta et all [20lll PLATO (|Rauer et al.l 120^. K 2 
(IHowell et al.ll2014n. SPECUL OOS (ICillon et ahlfeoTll 


SPIR ou (|Thibault et al.l[2012ll . and TESS (|Ricker et al.l 

[Ml). 


TABLE 9 

Candidates Accepted By Our Pipeline 


KID 

KOI 

Period (Days) 
Value - Err 

-1- Err 

Rp (R®) 
Value - Err 

-1- Err 

Value 

(F®) 

- Err 

-h Err Provenance^ 

2161536 

2130.01 

16.85595526 

6.07E-05 

5.96E-05 

1.93 

0.26 

0.30 

6.27 

1.56 

1.86 

0 

2556650 

2156.01 

2.85234932 

3.29E-06 

3.38E-06 

1.88 

0.22 

0.27 

37.94 

8.21 

9.58 

0 

2715135 

1024.01 

5.74773353 

8.53E-06 

8.58E-06 

1.49 

0.19 

0.21 

22.53 

5.97 

7.60 

0 

2973386 

3034.01 

31.02089195 

2.46E-04 

2.30E-04 

1.62 

0.33 

0.37 

2.31 

0.85 

1.12 

0 

3426367 

2662.01 

2.10434036 

4.05E-06 

4.03E-06 

0.56 

0.09 

0.09 

28.25 

8.02 

9.82 

0 

3642335 

3010.01 

60.86661711 

4.91E-04 

5.80E-04 

2.37 

0.27 

0.33 

0.94 

0.18 

0.21 

0 

3749365 

1176.01 

1.97376228 

4.23E-07 

4.24E-07 

9.11 

2.01 

2.01 

74.26 

36.13 

60.43 

0 

4061149 

1201.01 

2.75759481 

6.97E-06 

6.65E-06 

1.21 

0.18 

0.20 

39.76 

10.43 

12.72 

0 

4139816 

812.01 

3.34022436 

8.45E-06 

8.40E-06 

2.11 

0.32 

0.33 

38.49 

12.17 

16.10 

1 

4139816 

812.02 

20.05993080 

2.82E-04 

3.03E-04 

2.14 

0.33 

0.36 

3.53 

1.12 

1.48 

1 

4139816 

812.03 

46.18428829 

6.06E-04 

6.48E-04 

1.90 

0.29 

0.31 

1.16 

0.37 

0.49 

1 

4139816 

812.04 

7.82527952 

1.24E-04 

1.37E-04 

1.09 

0.18 

0.22 

12.37 

3.91 

5.19 

1 

4172805 

4427.01 

147.66174063 

1.51E-03 

2.16E-03 

1.56 

0.23 

0.25 

0.17 

0.05 

0.06 

0 

4832837 

605.01 

2.62811645 

2.03E-06 

2.04E-06 

2.52 

0.41 

0.43 

69.25 

28.58 

43.81 

0 

4832837 

605.02 

5.06549822 

4.22E-05 

4.37E-05 

0.97 

0.17 

0.19 

28.88 

11.92 

18.35 

0 

4913852 

818.01 

8.11438395 

1.25E-05 

1.24E-05 

2.14 

0.26 

0.27 

11.23 

2.81 

3.47 

0 

5364071 

248.01 

7.20387100 

8.00E-06 

8.00E-06 

2.25 

0.27 

0.27 

14.00 

3.64 

4.55 

3 

5364071 

248.02 

10.91273200 

2.10E-05 

2.10E-05 

2.74 

0.36 

0.45 

8.04 

2.08 

2.58 

3 

5364071 

248.03 

2.57654900 

3.00E-06 

3.00E-06 

1.55 

0.18 

0.18 

55.11 

14.21 

17.62 

3 

5364071 

248.04 

18.59610800 

7.90E-05 

7.90E-05 

1.31 

0.17 

0.17 

3.94 

1.01 

1.28 

3 

5384713 

3444.02 

60.32665084 

4.47E-05 

4.47E-05 

2.98 

0.34 

0.36 

1.08 

0.18 

0.21 

4 

5531953 

1681.00 

21.91384343 

1.83E-04 

1.93E-04 

1.03 

0.15 

0.17 

1.85 

0.47 

0.57 

0 


These distance estimates are based on the mean number of 
planets per star rather than the fraction of stars with planets. If 
potentially habitable planets are clustered such that M dwarfs host¬ 
ing potentially habitable planets typically feature more than one 


potentially habitable planet, then the distance estimates will need 
to be increased to account for the relatively flat nature of multi¬ 
planet systems orbiting M dwarfs dBallard fc Johnsonll201^b 
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TABLE 9 — Continued 


KID 

KOI 

Period (Days) 
Value - Err 

+ Err 

Rp (R®) 
Value - Err 

+ Err 

Fp (F®) 
Value - Err 

+ Err Provenance^ 

5531953 

1681.01 

6.93911381 

2.60E-05 

2.60E-05 

0.99 

0.13 

0.13 

8.55 

2.07 

2.38 

2 

5531953 

1681.02 

1.99281275 

7.47E-06 

7.21E-06 

0.72 

0.10 

0.11 

45.27 

11.46 

13.98 

0 

5531953 

1681.03 

3.53105829 

2.39E-05 

2.40E-05 

0.68 

0.10 

0.12 

21.12 

5.34 

6.52 

0 

5617854 

1588.01 

3.51749675 

3.96E-06 

3.94E-06 

1.21 

0.15 

0.16 

40.95 

9.81 

11.76 

0 

5640085 

448.01 

10.13962300 

2.20E-05 

2.20E-05 

1.77 

0.22 

0.23 

9.50 

2.35 

2.84 

3 

5640085 

448.02 

43.59579200 

1.25E-04 

1.25E-04 

2.48 

0.31 

0.32 

1.36 

0.33 

0.41 

3 

5794240 

254.01 

2.45524062 

l.OOE-07 

l.OOE-07 

11.00 

0.56 

0.56 

67.54 

9.34 

10.62 

3 

5809954 

1902.01 

137.86453830 

2.87E-04 

2.87E-04 

1.99 

0.40 

0.44 

0.23 

0.06 

0.08 

4 

6382217 

2036.01 

8.41102527 

2.59E-05 

2.61E-05 

1.51 

0.30 

0.31 

13.17 

4.79 

6.30 

0 

6382217 

2036.02 

5.79529807 

3.24E-05 

3.52E-05 

1.00 

0.20 

0.21 

21.68 

7.86 

10.38 

0 

6435936 

854.01 

56.05318853 

1.31E-03 

1.26E-03 

2.09 

0.29 

0.32 

0.65 

0.16 

0.20 

1 

6497146 

3284.01 

35.23301880 

2.24E-04 

2.27E-04 

1.01 

0.16 

0.17 

1.31 

0.38 

0.47 

0 

6666233 

2306.01 

0.51240811 

5.59E-07 

5.53E-07 

0.94 

0.13 

0.14 

541.12 

143.41 

173.03 

0 

6679295 

2862.01 

24.57535492 

1.37E-04 

1.37E-04 

1.60 

0.18 

0.20 

2.51 

0.56 

0.67 

0 

6773862 

1868.01 

17.76080479 

2.70E-05 

2.62E-05 

2.13 

0.21 

0.24 

5.69 

1.12 

1.31 

0 

6867155 

868.01 

235.99802060 

3.78E-04 

3.78E-04 

8.88 

1.18 

1.17 

0.14 

0.04 

0.05 

3 

7021681 

255.01 

27.52199799 

4.98E-05 

4.96E-05 

2.56 

0.31 

0.31 

2.27 

0.54 

0.65 

1 

7021681 

255.02 

13.60335797 

2.02E-04 

2.01E-04 

0.75 

0.11 

0.12 

5.79 

1.39 

1.65 

1 

7094486 

1907.01 

11.35011141 

2.30E-05 

2.33E-05 

1.98 

0.20 

0.22 

9.32 

1.90 

2.25 

0 

7135852 

875.01 

4.22097130 

3.02E-06 

3.04E-06 

2.58 

0.45 

0.46 

29.90 

18.13 

36.70 

0 

7287995 

877.01 

5.95489377 

8.57E-06 

8.63E-06 

2.06 

0.29 

0.30 

20.25 

6.52 

8.89 

1 

7287995 

877.02 

12.03993346 

2.85E-05 

2.92E-05 

1.87 

0.26 

0.26 

7.92 

2.56 

3.48 

1 

7287995 

877.03 

20.83776773 

1.97E-04 

1.94E-04 

1.05 

0.16 

0.19 

3.81 

1.23 

1.67 

1 

7304449 

1702.01 

1.53818130 

3.09E-06 

3.12E-06 

0.84 

0.17 

0.18 

27.77 

10.07 

13.20 

0 

7447200 

676.01 

7.97251347 

1.82E-06 

1.82E-06 

2.88 

0.33 

0.32 

14.41 

3.22 

3.70 

6 

7447200 

676.02 

2.45323590 

4.75E-07 

4.75E-07 

3.67 

0.41 

0.42 

69.37 

15.43 

18.00 

6 

7455287 

886.01 

8.01026000 

3.00E-05 

3.00E-05 

2.39 

0.35 

0.43 

10.04 

2.56 

3.09 

3 

7455287 

886.02 

12.07238900 

l.OOE-04 

l.OOE-04 

1.28 

0.17 

0.17 

5.81 

1.49 

1.78 

3 

7455287 

886.03 

20.99569400 

1.43E-04 

1.43E-04 

1.40 

0.18 

0.19 

2.77 

0.70 

0.85 

3 

7603200 

314.01 

13.78113900 

1.10E-05 

l.lOE-05 

1.34 

0.16 

0.16 

5.98 

1.42 

1.68 

3 

7603200 

314.02 

23.08871300 

3.10E-05 

3.10E-05 

1.30 

0.16 

0.18 

3.01 

0.71 

0.84 

3 

7603200 

314.03 

10.31236400 

3.60E-05 

3.60E-05 

0.59 

0.07 

0.08 

8.81 

2.10 

2.47 

3 

7870390 

898.01 

9.77042372 

3.02E-05 

3.20E-05 

2.38 

0.29 

0.31 

10.52 

2.56 

3.13 

1 

7870390 

898.02 

5.16981333 

2.29E-05 

2.19E-05 

1.78 

0.22 

0.24 

24.57 

5.99 

7.28 

1 

7870390 

898.03 

20.09010000 

1.33E-04 

1.40E-04 

2.04 

0.25 

0.27 

4.02 

0.98 

1.20 

1 

7871954 

1515.01 

1.93703537 

2.77E-06 

2.71E-06 

0.95 

0.12 

0.13 

94.54 

22.56 

26.94 

0 

7871954 

1515.02 

7.06117534 

1.34E-05 

1.38E-05 

1.16 

0.15 

0.17 

16.86 

4.04 

4.81 

0 

7907423 

899.01 

7.11369666 

2.66E-05 

2.68E-05 

1.25 

0.18 

0.20 

8.49 

2.37 

2.90 

1 

7907423 

899.02 

3.30656751 

1.39E-05 

1.42E-05 

0.99 

0.15 

0.16 

23.55 

6.58 

8.03 

1 

7907423 

899.03 

15.36834908 

1.08E-04 

1.12E-04 

1.17 

0.17 

0.18 

3.04 

0.85 

1.04 

1 

8013419 

901.01 

12.73263426 

5.97E-06 

6.01E-06 

5.02 

0.86 

0.87 

7.98 

5.02 

10.44 

0 

8018547 

902.01 

83.94017967 

1.49E-03 

1.38E-03 

5.23 

0.76 

0.80 

0.62 

0.18 

0.22 

1 

8120608 

571.01 

7.26737224 

1.93E-05 

1.98E-05 

1.26 

0.17 

0.17 

11.52 

3.00 

3.63 

1 

8120608 

571.02 

13.34295116 

4.25E-05 

4.45E-05 

1.28 

0.17 

0.18 

5.12 

1.34 

1.61 

1 

8120608 

571.03 

3.88677934 

9.86E-06 

9.88E-06 

1.06 

0.14 

0.15 

26.51 

6.90 

8.35 

1 

8120608 

571.04 

22.40788397 

1.20E-04 

1.27E-04 

1.18 

0.16 

0.18 

2.57 

0.67 

0.81 

1 

8120608 

571.05 

129.94188177 

2.32E-03 

2.34E-03 

1.02 

0.14 

0.16 

0.25 

0.06 

0.08 

0 

8167996 

1867.01 

2.54954597 

1.25E-05 

5.91E-06 

1.12 

0.12 

0.14 

53.87 

11.96 

14.43 

1 

8167996 

1867.02 

13.96935372 

9.90E-05 

1.07E-04 

1.81 

0.21 

0.26 

5.58 

1.24 

1.49 

1 

8167996 

1867.03 

5.21231974 

5.15E-05 

5.22E-05 

1.16 

0.13 

0.15 

20.77 

4.61 

5.56 

1 

8189801 

2480.01 

0.66682583 

9.30E-07 

9.58E-07 

1.31 

0.23 

0.25 

466.69 

145.16 

180.18 

0 

8229458 

2238.01 

1.64680052 

2.72E-06 

2.81E-06 

0.93 

0.10 

0.12 

118.41 

22.88 

26.30 

0 

8235924 

2347.01 

0.58800191 

7.34E-07 

7.18E-07 

1.02 

0.10 

0.13 

550.07 

101.56 

116.51 

0 

8351704 

1146.01 

7.09712741 

2.98E-05 

3.27E-05 

0.82 

0.13 

0.14 

7.97 

2.35 

2.94 

0 

8367644 

1879.01 

22.08557431 

5.05E-05 

4.98E-05 

2.42 

0.37 

0.39 

1.96 

0.57 

0.71 

0 

8509442 

2992.01 

82.65995306 

6.70E-04 

7.23E-04 

2.23 

0.35 

0.37 

0.70 

0.20 

0.25 

0 

8547140 

1266.01 

11.41929554 

2.42E-05 

2.45E-05 

1.54 

0.26 

0.27 

7.27 

2.55 

3.52 

0 

8561063 

961.01 

1.21377022 

3.76E-07 

3.76E-07 

0.87 

0.21 

0.21 

17.83 

7.89 

11.04 

2 

8561063 

961.02 

0.45328734 

8.60E-08 

8.60E-08 

0.90 

0.21 

0.21 

67.88 

31.62 

54.63 

3 

8561063 

961.03 

1.86511477 

1.12E-06 

1.12E-06 

0.69 

0.16 

0.16 

10.00 

4.38 

6.24 

2 

8631751 

2453.01 

1.53051487 

1.96E-06 

1.88E-06 

1.13 

0.21 

0.23 

61.55 

21.28 

28.63 

0 

8845205 

463.01 

18.47717612 

1.92E-04 

2.26E-04 

1.67 

0.28 

0.28 

1.27 

0.41 

0.51 

1 

8874090 

1404.01 

13.32391953 

5.44E-05 

5.62E-05 

1.35 

0.27 

0.28 

4.79 

2.11 

3.29 

0 

8874090 

1404.02 

18.90623146 

2.81E-04 

3.04E-04 

0.96 

0.20 

0.21 

3.01 

1.32 

2.06 

0 

8890150 

2650.01 

34.98870981 

2.25E-03 

3.62E-03 

0.96 

0.15 

0.17 

1.17 

0.38 

0.62 

1 

9214942 

1403.01 

18.75471931 

5.57E-05 

5.57E-05 

1.78 

0.28 

0.29 

5.03 

1.59 

2.11 

0 

9388479 

936.01 

9.46787010 

6.62E-05 

6.29E-05 

2.18 

0.30 

0.31 

6.26 

1.68 

2.04 

1 

9388479 

936.02 

0.89303356 

4.90E-06 

4.81E-06 

1.26 

0.18 

0.18 

145.89 

39.21 

47.53 

1 

9390653 

249.01 

9.54930070 

1.50E-05 

1.42E-05 

1.62 

0.22 

0.22 

4.54 

1.20 

1.44 

1 

9427402 

1397.01 

6.24703152 

1.02E-05 

1.03E-05 

2.23 

0.23 

0.26 

21.87 

4.24 

4.89 

0 

9573685 

2057.01 

5.94565639 

1.53E-05 

1.52E-05 

1.19 

0.17 

0.19 

21.39 

5.60 

6.80 

0 

9575728 

5692.01 

2.64180062 

1.99E-05 

2.22E-05 

0.46 

0.05 

0.06 

64.99 

13.21 

15.76 

0 

9710326 

947.01 

28.59885823 

2.34E-04 

2.18E-04 

1.92 

0.26 

0.27 

1.81 

0.47 

0.57 

1 

9757613 

250.01 

12.28293091 

1.43E-05 

1.41E-05 

2.89 

0.44 

0.45 

7.83 

2.76 

3.88 

1 

9757613 

250.02 

17.25116937 

3.21E-05 

3.18E-05 

2.51 

0.40 

0.45 

4.98 

1.75 

2.46 

1 

9757613 

250.03 

3.54391419 

1.15E-05 

1.15E-05 

1.06 

0.17 

0.18 

41.08 

14.46 

20.34 

0 






23 


The Occurrence Rate of Potentially Habitable Planets Orbiting M dwarfs 

TABLE 9 — Continued 


KID 

KOI 

Period (Days) 
Value - Err 

-1- Err 

Rp (R®) 
Value - Err 

-1- Err 

} 

Value 

"’p (F®) 

- Err 

-h Err Provenance^ 

9757613 

250.04 

46.82763492 

1.60E-04 

1.62E-04 

2.10 

0.33 

0.36 

1.31 

0.46 

0.65 

1 

9787239 

952.01 

5.90129534 

2.03E-05 

1.94E-05 

2.13 

0.27 

0.29 

16.39 

4.14 

5.11 

1 

9787239 

952.02 

8.75208025 

2.04E-05 

2.09E-05 

1.98 

0.27 

0.34 

9.70 

2.46 

3.02 

0 

9787239 

952.03 

22.78068566 

1.28E-04 

1.31E-04 

2.30 

0.29 

0.29 

2.71 

0.69 

0.84 

1 

9787239 

952.04 

2.89601166 

8.92E-06 

9.34E-06 

1.09 

0.14 

0.14 

42.36 

10.69 

13.21 

0 

9787239 

952.05 

0.74295972 

3.06E-06 

3.65E-06 

0.91 

0.12 

0.13 

260.04 

66.04 

80.88 

1 

10027247 

2418.01 

86.82818564 

1.57E-03 

1.27E-03 

1.24 

0.17 

0.23 

0.35 

0.09 

0.11 

0 

10027323 

1596.01 

5.92367766 

1.81E-05 

1.89E-05 

1.04 

0.16 

0.16 

18.61 

5.52 

7.08 

0 

10027323 

1596.02 

105.35789935 

5.05E-04 

5.36E-04 

1.90 

0.28 

0.30 

0.40 

0.12 

0.15 

0 

10073672 

2764.01 

2.25297009 

7.20E-06 

7.26E-06 

1.30 

0.22 

0.24 

85.43 

27.60 

35.63 

0 

10118816 

1085.01 

7.71790094 

5.17E-05 

4.78E-05 

0.94 

0.14 

0.16 

15.00 

4.82 

6.54 

0 

10166274 

1078.01 

3.35372553 

1.04E-05 

1.06E-05 

1.86 

0.25 

0.26 

34.92 

9.20 

11.35 

1 

10166274 

1078.02 

6.87749137 

2.35E-05 

2.40E-05 

2.10 

0.28 

0.31 

13.39 

3.51 

4.35 

1 

10166274 

1078.03 

28.46416009 

2.41E-04 

2.27E-04 

1.93 

0.26 

0.29 

2.02 

0.53 

0.65 

1 

10329835 

2058.01 

1.52372638 

2.40E-06 

2.46E-06 

1.05 

0.14 

0.16 

131.36 

34.28 

41.60 

0 

10332883 

1880.01 

1.15116708 

6.89E-07 

6.97E-07 

1.30 

0.20 

0.21 

181.87 

53.73 

66.69 

0 

10340423 

736.01 

18.79420980 

5.75E-05 

5.95E-05 

1.82 

0.30 

0.31 

2.81 

0.89 

1.14 

0 

10340423 

736.02 

6.73899074 

3.10E-05 

3.13E-05 

1.17 

0.19 

0.20 

11.03 

3.48 

4.48 

0 

10386984 

739.01 

1.28707954 

4.42E-06 

4.48E-06 

1.45 

0.18 

0.19 

129.24 

30.84 

36.80 

1 

10388286 

596.01 

1.68269527 

5.88E-06 

5.95E-06 

1.38 

0.19 

0.22 

76.32 

19.27 

23.05 

1 

10489206 

251.01 

4.16438060 

2.78E-06 

2.77E-06 

2.61 

0.31 

0.31 

30.10 

6.98 

8.25 

1 

10489206 

251.02 

5.77446840 

6.66E-05 

6.90E-05 

0.90 

0.12 

0.15 

19.47 

4.52 

5.33 

1 

10525027 

2006.01 

3.27346499 

5.99E-06 

6.07E-06 

0.74 

0.10 

0.11 

35.43 

9.80 

12.56 

0 

10525049 

4252.01 

15.57135826 

9.44E-05 

9.36E-05 

0.73 

0.14 

0.15 

5.54 

2.00 

2.61 

0 

10591855 

2845.01 

1.57408931 

5.05E-06 

5.08E-06 

0.89 

0.10 

0.13 

135.82 

26.43 

30.77 

0 

10670119 

2179.01 

14.87155876 

5.40E-05 

5.01E-05 

1.27 

0.19 

0.20 

3.15 

0.90 

1.10 

0 

10670119 

2179.02 

2.73277006 

5.46E-06 

5.49E-06 

1.08 

0.16 

0.18 

30.19 

8.59 

10.50 

0 

10717241 

430.01 

12.37646610 

1.30E-05 

1.30E-05 

2.04 

0.26 

0.27 

6.15 

1.57 

1.92 

3 

10717241 

430.02 

9.34052690 

8.42E-05 

8.42E-05 

0.73 

0.10 

0.10 

8.95 

2.30 

2.78 

3 

10925104 

156.01 

8.04133834 

9.01E-06 

8.98E-06 

1.39 

0.18 

0.20 

15.72 

3.58 

4.24 

1 

10925104 

156.02 

5.18856137 

9.66E-06 

9.67E-06 

1.04 

0.13 

0.15 

28.20 

6.40 

7.60 

1 

10925104 

156.03 

11.77614238 

7.55E-06 

7.52E-06 

2.02 

0.23 

0.23 

9.45 

2.15 

2.55 

0 

11129738 

1427.01 

2.61301749 

5.76E-06 

5.70E-06 

1.36 

0.20 

0.22 

58.27 

17.40 

22.67 

0 

11187837 

252.01 

17.60462615 

3.24E-05 

3.21E-05 

2.36 

0.29 

0.29 

3.79 

0.93 

1.10 

1 

11192235 

2329.01 

1.61535973 

2.53E-06 

2.53E-06 

1.17 

0.22 

0.24 

101.54 

20.26 

23.51 

2 

11348997 

2090.01 

5.13248495 

9.13E-06 

9.09E-06 

1.52 

0.25 

0.26 

18.98 

5.92 

7.38 

0 

11497958 

1422.01 

5.84164124 

1.02E-05 

1.02E-05 

3.49 

0.30 

0.50 

14.89 

5.58 

8.77 

5 

11497958 

1422.02 

19.85028393 

5.85E-05 

5.85E-05 

3.10 

0.51 

1.11 

2.94 

1.30 

1.78 

5 

11497958 

1422.03 

10.86443187 

4.97E-05 

4.97E-05 

2.34 

0.34 

0.87 

5.09 

1.96 

7.23 

5 

11497958 

1422.04 

63.33666060 

5.57E-04 

5.57E-04 

2.30 

0.34 

0.66 

0.61 

0.37 

0.53 

5 

11497958 

1422.05 

34.14189000 

2.64E-04 

2.64E-04 

2.12 

0.34 

0.79 

1.42 

0.72 

0.93 

5 

11752906 

253.01 

6.38316009 

9.25E-06 

9.25E-06 

2.68 

0.39 

0.40 

17.61 

6.43 

9.30 

1 

11752906 

253.02 

20.61727169 

2.66E-04 

2.85E-04 

1.48 

0.23 

0.25 

3.69 

1.35 

1.94 

1 

11754553 

775.01 

16.38481298 

7.05E-05 

7.02E-05 

1.90 

0.26 

0.27 

5.56 

1.61 

2.08 

1 

11754553 

775.02 

7.87740709 

1.12E-05 

1.12E-05 

2.10 

0.28 

0.31 

14.75 

4.23 

5.53 

3 

11754553 

775.03 

36.44516433 

2.60E-04 

2.91E-04 

1.96 

0.27 

0.29 

1.91 

0.55 

0.72 

1 

11768142 

2626.01 

38.09723780 

2.86E-04 

2.86E-04 

2.36 

0.53 

0.44 

0.92 

0.46 

0.57 

5 

11852982 

247.01 

13.81496561 

6.45E-05 

6.47E-05 

1.61 

0.21 

0.22 

5.22 

1.27 

1.50 

1 

11853130 

3263.01 

76.87935073 

4.79E-05 

4.79E-05 

6.83 

1.18 

1.27 

0.31 

0.08 

0.09 

2 

11853255 

778.01 

2.24336795 

1.13E-05 

l.llE-05 

1.37 

0.19 

0.20 

52.26 

14.75 

18.85 

1 

11923270 

781.01 

11.59822478 

1.50E-05 

1.49E-05 

2.82 

0.37 

0.39 

6.10 

1.55 

1.90 

0 

12066335 

784.01 

19.27103179 

3.52E-04 

3.51E-04 

1.81 

0.25 

0.28 

3.29 

0.88 

1.10 

1 

12066335 

784.02 

10.06525147 

2.48E-05 

2.46E-05 

1.54 

0.20 

0.22 

7.83 

2.10 

2.63 

0 

12066569 

3282.01 

49.27623343 

4.37E-04 

5.12E-04 

2.07 

0.29 

0.31 

1.30 

0.34 

0.41 

0 

12302530 

438.01 

5.93119294 

5.02E-06 

5.02E-06 

1.80 

0.24 

0.25 

23.47 

7.53 

10.42 

3 

12302530 

438.02 

52.66160633 

1.53E-04 

1.53E-04 

1.80 

0.25 

0.27 

1.28 

0.41 

0.56 

0 

12352520 

3094.01 

4.57700369 

1.33E-05 

1.34E-05 

1.39 

0.17 

0.20 

24.39 

5.04 

5.86 

0 

12506770 

1577.01 

2.80624351 

8.14E-06 

8.33E-06 

1.40 

0.23 

0.25 

59.53 

18.67 

24.43 

0 


^ The planet parameters are from our fits to the long cadence 
data (provenance = 0), our fits to the short cadence data (prove¬ 
nance = 1), the NASA Exoplanet Archive (provenance = 2), 
IRowe et alJ 1)20141 . provenance = Sh lSwift et al.l(|2015l . provenance 
= 4), ICartier et al.l (120141 . provenance = 5), and lloannidis et al.l 
lf20Tl . provenance = 6). 
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